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ABSTRACT 

We present semi-analytic models of galactic outflows that are constrained by available ob- 
servations on high redshift star formation and reionization. Galactic outflows are modeled in 
a manner akin to models of stellar wind blown bubbles. Large scale outflows can generically 
escape from low mass halos {Ad < 10^ Mq) for a wide range of model parameters while 
this is not the case in high mass halos (M > 10^^ Mq). The flow generically accelerates 
within the halo virial radius, then starts to decelerate, and traverses well into the intergalactic 
medium (IGM), before freezing to the Hubble flow. The acceleration phase can result in shell 
fragmentation due to the Rayleigh-Taylor instability, although the final outflow radius is not 
significantly altered. The gas phase metallicity of the outflow and within the galaxy are com- 
puted assuming uniform instantaneous mixing. Ionization states of different metal species are 
calculated and used to examine the detectability of metal fines from the outflows. The global 
influence of galactic outflows is also investigated using porosity weighted averages and prob- 
ability density functions of various physical quantities. Models with only atomic cooled halos 
significantly fifl the IGM at z 3 with metals (with -2.5 > [Z/Zq] > -3.7), the actual 
extent depending on the efficiency of winds, the initial mass fimction (IMF) and the frac- 
tional mass that goes through star formation. The reionization history has a significant effect 
on the volume filling factor, due to radiative feedback. In these models, a large fraction of 
outflows at z ^ 3 are supersonic, hot (T > lO^K) and have low density, making metal 
lines difficult to detect. They may also result in significant perturbations in the IGM gas on 
scales probed by the Lyman-a forest. On the contrary, models including molecular cooled 
halos with a normal mode of star formation can potentially volume fill the universe at z > 8 
without drastic dynamic effects on the IGM, thereby setting up a possible metallicity floor 
(—4.0 < [Z/Zq] < —3.6). The order unity fluctuations at z ^ 8 that becomes the mildly 
non-linear fluctuations traced by Lyman-a forest at z < 4 will then have this metaUicity. In- 
terestingly, molecular cooled halos with a "top-heavy" mode of star formation are not very 
successful in establishing the metalUcity floor because of the additional radiative feedback, 
that they induce. 

Key words: cosmology: theory - early Universe - galaxies : formation - evolution - high- 
redshift - intergalactic medium - stars: winds, outflows 



1 INTRODUCTION 

The rapid growth of observations of the high redshift universe has 
raised several intriguing questions regarding physics of galaxy for- 
mation and the physical state of the intergalactic medium (IGM). 
Some of the important issues are: how and when the dark ages 
ended with the reionization of the IGM, the origin of the metals, 
and temperature of the low density IGM traced by Lyman-a forest 
at z ~ 2.5. In this work we concentrate on the issue of IGM metal 
enrichment. 
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Presence of metals in the Lyman-a forest (with r > 1) is 
now well established though C iv (Tytler et al. 1995; Songaila 
& Cowie, 1996) and O vi (Carswell et al. 2002; Simcoe et al. 
2002; Bergeron et al. 2002) absorption lines detected in the echelle 
specU-a of high redshift QSOs. The measured N(C iv)/N(H I) at 
a ~ 3 are consistent with [C/H]~ —2.5 with large errors (see 
Rauch et al. 1997). Songaila (2001) has shown that the C IV col- 
umn density distribution function is consistent with being invari- 
ant between z = 1.5 and z — 5 and a minimum metallicity of 
[C/H]~ —3.5 is already in place at z ~ 5.5 (see also Ryan- Weber 
et al. 2006). As under-dense space occupy most of the volume, 
measuring metallicity in the low density region is very important. 
Given the expected low metallicities and large ionization correc- 
tions, direct detection of metals from such low neutral hydrogen 
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optical depth (i.e 13 < log [N(H I) cm ] < 14) is currently im- 
possible and usually the metallicity is estimated using pixel statis- 
tics ( [Ellison et al. 2000[ [Schaye et al. 2003[ [Aracil et al. 2004| l. 
Schaye et al. (2003) obtained [C/H]~ -3.5 for log 5 > -0.5 
at 2 = 1.8 — 4.1. Here, 5 being the over density defined as 
5 = p/p > 10{N{H i)/10i5cm-2)2/3[(l + z)/4]-^ (Schaye 
2001) and p is the mean IGM density. However, Aracil et al. (2004) 
have not detected C IV from r < 1 Lyman-a systems in their 
UVES data. Thus, what fraction of the IGM is filled with metals is 
a subject of ongoing debates and investigations (Schaye et al. 2003; 
Aracil et al. 2004; Scannapieco et al. 2006). In any case, in the stan- 
dard framework of LCDM models the Lyman-a forest absorption 
lines originate from density fluctuations that are either in the linear 
or in the quasi-linear regime (Bi & Davidsen, 1997). Thus, metals 
that may be present in these regions need to be transported from 
the neighboring star-forming regions. The amount and distribution 
of metals in the Lyman-a forest provides information on different 
feedbacks from star forming galaxies. 

Observations of high-z Lyman break galaxies (LBGs) fre- 
quently show galactic scale superwinds (Pettini et al. 2001), high 
value for the escape fraction of the UV photons (Steidel et al. 2001) 
and a strong correlation between the C IV absorption systems and 
the LBGs at 2 ^ 3 within an impact parameter of ~ 50 kpc (Adel- 
berger et al. 2005). Profiles of C IV and O VI absorption lines seen 
in the high redshift damped Lyman-a systems (DLAs) are also con- 
sistent with them originating from outflows in DLA protogalaxies 
(Fox et al. 2007, 2007a). Adelberger et al. (2005) have argued that 
considerable fraction of high column density C IV systems may 
originate from the large scale galactic outflows. Clustering proper- 
ties of strong C IV absorption lines (Scannapieco et al. 2006) and 
the velocity profiles C IV absorption lines (Songaila et al. 2006) are 
consistent with a good fraction of the high column density C IV sys- 
tems originating from the region near massive galaxies. However, 
the important question is whether these metals are freshly emitted 
by the bright galaxies or due to biased clustering of gas ejected 
by low mass galaxies from early (say 2 > 6) epochs (Porciani & 
Madau, 2005). Aguirre et al. (2005) by comparing spectra predicted 
by various simulations that includes winds found that the predicted 
metal distribution in the models are highly inhomogeneous and can 
not reproduce the observed probability distribution of C IV absorp- 
tion. They suggested that strong winds from galaxies at 2 < 6 
cannot fully explain the observed enrichment and additional pre- 
enrichment from higher redshift galaxies are needed. 

Note that star formation is the key element in controlling the 
outflows from galaxies. Star formation also has an important ef- 
fect on the reionization of the universe which in turn affects the 
star formation in low mass halos through radiative feedback. There 
exists a growing body of data constraining the star formation rate 
(SFR) density in the high redshift universe (jBou wens et al. 2005[ 
[Richards et al. 2 006 , Hopkins & Beacom 2006 1. Also observational 
constraints on the reionization are provided by the spectra of high- 
est redshift QSOs (Fan et al. 2006) and ongoing WMAP satellite 
observations of Cosmic Microwave Background (CMB) polariza- 
tion (Spergel et al. 2007). Therefore one must deal with reionization 
and galactic outflows in a manner consistent with the observed star 
formation history of the universe. In our previous paper Samui, Sri- 
anand & Subramanian (2007; here after Paper I) we have described 
a set of models which correctly produce the observed luminosity 
functions of LBGs and hence the SFR density at different redshifts. 
Here we take those constrained models of star formation and reion- 
ization from Paper I and used them to predict the various properties 
galactic outflows and their impact on the IGM such as volume fill- 



ing factor of wind/metals, mean metallicity, temperature of these 
polluted region etc. 

There are several attempts to model the galactic outflow us- 
ing semi-analytic calculations (for example, Madau, Ferrara & 
Rees 2001; Scannapieco, Ferrara & Madau 2002; Furlanetto & 
Loeb 2003; Scannapieco 2005). The model presented here will be 
broadly in line with these above studies with significant improve- 
ments. In particular, the star formation is continuous and duration 
of the star formation and the amount of gas going into stars are 
constrained by the observed high-2 luminosity function (see Paper 
I and section 2 below). This is very different from the previous at- 
tempts where the star formation is usually in the form of bursts. Our 
model takes into account the reionization and the radiative feedback 
to the star formation in a self-consistent way. This is very important, 
in particular, to estimate the metal pollution due to low mass halos 
and constrain the initial mass function (IMF) at high redshifts. We 
discuss various other improvements and their effects in detail. 

The paper is organized in following manner. In the next sec- 
tion we first briefly outline the model of star formation. In sec- 
tion 3 we introduce our model of galactic outflows in detail. The 
structural properties of the outflow and sensitivity of our results to 
the adopted initial conditions are discussed in detail in section 4. 
In section 5 we discuss the dependence of outflow properties on 
the model parameters. In section 6 we investigate the growth of 
Rayleigh-Tayler instabilities in the accelerating phase of the out- 
flow and study its consequences. The ionization correction and the 
detectability of the winds in different stages of its evolution are dis- 
cussed in section 7. The global effects of the outflow are studied in 
section 8 and the summary and conclusions are provided in section 
9. In most of this work we use the cosmological parameters consis- 
tent with the recent WMAP data (Q = 1, Sl„ = 0.26, Qa = 0.74, 
fib = 0.044, h = 0.71, as = 0.75 and Us = 0.95). 



2 STAR FORMATION RATE 

We follow the same prescription as in Samui, Srianand & Subrama- 
nian (2007) to get the star formation rate (SFR) density in a semi- 
analytical fashion. Here, we only briefly outline the model. We use 
the modified Press-Schechter (PS) formalism of Sasaki (1994) to 
calculate the number density of collapsed objects having mass be- 
tween M and M + dM, which are formed at the redshift interval 
{zc, Zc + dzc) and survive till redshift 2. This is given by (Chiu & 
Ostriker 2000; Choudhury & Srianand 2002) 



N{M,z,Zc)dM dzc = Nm{z^) 
D{zc) 



D{zc 

dZr. 



D{zc 



D(Z) H(2c)(l+2,) 



dM. 



(1) 



Here, Nm{zc) dM is the PS mass function 1 Press & Schechter 
[I974[ l, 5c is the critical over density for collapse, usually taken to 
be equal to 1.686. Further, H{z) is the Hubble parameter, D{z) 
the growth factor for linear perturbations and cr(A/) the rms mass 
fluctuation at a mass scale M. 

We model the star formation in a given halo of mass M col- 
lapsed at Zc and observed at 2 as. 



MsAM,z,Zc) = /* 



t{z) - t{Zc 



X exp 



dyn 
t{z) - t{Zc) 



(2) 
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Here, /* is the fraction of total baryonic mass in a halo that will be 
converted to stars. The function t{z) gives the age of the universe 
at redshift z; thus, t{z) — t{zc) is the age of the collapsed halo 
at z and tdyn is the dynamical time-scale ^Chiu & Ostriker 2000| 
[Barkana & Loeb 2001) . As mentioned in Paper I the duration of 
star formation activity in a halo depends on the value of k. Note 
that our prescription of SFR as a function of time (i.e. Eq. |2j is 
purely empirical and not obtained by taking into account all possi- 
ble physical processes in the interstellar medium (ISM). However, 
as the model parameters (/, & n) are constrained by the observed 
luminosity functions, we can treat it as the net SFR as a function 
of time resulting from various competing physical processes that 
govern the star formation. 

Having modeled the formation rate of halos and the star for- 
mation in an individual halo, we can calculate the comoving star 
formation rate density as, 

oo oo 

Psf{z)^ J dzc J dM'MsF{M',z,Zc)xN{M',z,Zc). (3) 

The lower mass cutoff (Miow) at a given epoch is decided by the 
cooling efficiency of the gas and different feedback processes (see 
Paper I for detail discussions). We consider models with Afiow cor- 
responding to a virial temperature, T^ir ~ 10* K ( as "atomic cool- 
ing model") and 300 K ("molecular cooling model") for the neutral 
gas. For ionized regions of the universe, our models assume com- 
plete suppression of star formation in halos below circular velocity 
Vc ~ 35 km s^^, no suppression above circular velocity of 95 km 
s~^ and a linear fit from 1 to for the intermediate masses (as in 
Bromm & Loeb 2002). The SFR in the high mass halos are re- 
duced by a suppression factor [1 + {M /W^'^ Mq)^]~^ . However, 
our models at present do not incorporate the chemical and SNe 
feedback. 

The reionization history of the universe is also calculated as 
outlined in Paper I. We assume all the Lyman-continuum pho- 
ton which escape the star forming region (with escape fraction 
fesc = 0.1) are used in ionization. We only consider the ionization 
of hydrogen and take case B recombination (with a recombination 
coefficient as) to calculate the recombination rate. The redshift 
evolution of ionized hydrogen fraction {fun) is given by. 



A^^ At 



Az 



nuiz) dz 



aBnH{z)fHiiC 



dt 
di' 



(4) 



Here, is the rate of UV photons escaping into the IGM and 
obtained from psF (z). Further, hh (z) is the proper number density 
of the hydrogen atoms and C is the clumping factor of the IGM. 
The redshift evolution of C is assumed to be in the form C{z) — 
1 + 9 X [7/(1 + z)f for z > 6 and C = 10 for z < 6. (Haiman 
& Bryan 2006) A'^^ for a given IMF is calculated as explained in 
Paper I (see Eq. (13) and Table 1 there). 

In Paper I, the observed UV luminosity function of galaxies at 
3 < 2 < 6 are well fitted by our models with /, =0.5 and k = 1 
for the cosmological parameters constrained by WMAP 3rd year 
data. Models with /, less than 0.5 require k < 1. However, the 
measured ages of high-z star formation galaxies favor k ~ 1. It has 
also been pointed out in Paper I, that models with slightly higher 
values of as and ris than that derived from the 3rd year WMAP 
data require lower values of /, even for k = 1. Thus, in our global 
models discussed in the following sections we will use three sets of 
parameters for the star formation models (i) k = 1 and /, = 0.5; 

(ii) K = 0.5 and /. = 0.25 with as = 0.75 & = 0.95 and 

(iii) K = 1 and /* = 0.25 with as = 0.85 and Us = 1.00. All 



these models reproduce the observed high-z luminosity functions 
reasonably well. Models with low values of k and /* correspond to 
cases where the fraction of cold gas that can be used for forming 
stars decreases more rapidly as a function of time. 

In Paper I, it has been shown that various ongoing deep imag- 
ing surveys will not be able to directly detect the molecular cooled 
halos that may be present prior to reionization. However, presence 
of these objects influences the reionization history and hence the 
extent of radiative feedback at jz > 6 (see section 6 in Paper I). 
As the gravitational potential of these sources are expected to be 
small it is usually believed that even a small amount of star forma- 
tion in these objects can drive outflows. Therefore, presence of star 
formation in molecular cooled halos can have an important role to 
play in the pre-enrichment of the IGM to some metallicity floor We 
discuss this issue by considering the models with molecular cooled 
halos that consistently produce reionization history constrained by 
the WMAP 3rd year data (see the list in Table 3 of Paper I). 



3 OUTFLOWS FROM GALAXIES INTO THE IGM 

The metals detected in the IGM can only have been synthesized by 
stars in galaxies, and galactic outflows are the primary means by 
which they can be transported from galaxies into the IGM. The me- 
chanical energy that drives such outflows may arise either from an 
active galactic nuclei (AGN) in the galaxy or from the supemovae 
(SNe) explosions associated with the star formation activities in the 
galaxy. Here we concentrate on the effects of the star formation ac- 
tivity and the resulting SNe in high redshift protogalaxies. 



3.1 The general outflow scenario 

When a single SNe explodes, it creates a bubble of shock heated 
interstellar medium (ISM) around itself which expands superson- 
ically into surrounding ISM. The occurrence of clustered and co- 
herent explosions of several SNe can lead to the merger of the as- 
sociated supernova remnants to form a super bubble. This super 
bubble expands as it is fed by consequent SNe explosions, and de- 
celerates as it sweeps up the ISM of the galaxy. In the case of a 
disk galaxy, as the super bubble radius approaches the disk scale 
height, the decreasing density of the ISM generally causes a reac- 
celeration of the swept up ISM shell, which then begins to fragment 
through growing Rayleigh-Taylor (RT) instabilities (see the discus- 
sion in Veilleux, Cecil & Bland-Hawthorn 2005). These fragments 
and the shock heated hot gas are then vented out into the galaxy 
halo. Further, in a protogalaxy with high enough rates of star for- 
mation, supernova remnants, even if randomly distributed, could 
fill a significant volume of the whole galaxy to create a galaxy wide 
super bubble. This will also feed hot gas into the galactic halo. 

In the absence of any gas in the halo of the galaxy, the hot 
gas being fed into the halo by super bubbles would escape as a 
thermally driven wind into the intergalactic space. However, it is 
more likely that the halo of a forming galaxy itself has residual gas 
which has not yet collapsed to the centre or gas which is continuing 
to fall in from the intergalactic medium. For example, if accretion 
of gas is bimodal, with both an early dominant cold mode along 
filaments and a latter hot mode via an accretion shock (cf. Dekel & 
Birnboim 2006; Keres et al. 2005), then significant star formation 
can occur before the hot gas in the halo completely accretes onto 
the galaxy. In this case, the bubble of SNe heated hot gas ejected 
from the galaxy, will be initially confined by the outside medium 
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Figure 1. Schematic diagram of tiie structure of tlie galactic outflow in our 
models. Panel (a) is for the 'pressure driven flow' and panel (b) is for the 
'momentum driven flow'. 



the fresh supply of SNe heated gas from the galaxy ■would be at 
a higher pressure. A sufficiently large pressure difference between 
the SNe heated hot gas and the bubble gas would lead to a thermally 
driven wind being driven out of the galaxy, where the thermal en- 
ergy of the SNe heated gas is converted into directed kinetic energy 
of the wind (cf. Chevalier & Clegg 1985). Such a wind will be again 
stopped at a distance where its ram pressure equals the previously 
created hot bubble pressure. This is the origin of the inner shock at 
radius i?i . The above evolutionary sequence is different from the 
case of a stellar wind blown bubble, where the wind from the star 
has already a high velocity and low temperature that it leads to the 
"reverse shock" at i?i right from the early stages. Nevertheless the 
subsequent evolution can be very similar to the stellar case. 

In the case when the hot bubble material in region (b) between 
Ri and Rc cools efficiently, Ri can approach Rc with zone (b) be- 
coming of negligible thickness. The wind material from the galaxy 
then directly impinges on the shell (as shown in panel (b) in Fig.[T|l. 
The wind momentum is directly transferred to the shell and the out- 
flow then becomes momentum driven. We shall refer to the outflow 
then as a "momentum driven outflow". 

Note, the thin shell approximation we will be using is valid if 
the cooling time in post-shocked halo/IGM gas is small compared 
to the expansion time of the bubble. For interstellar bubbles, be- 
ing blown into a constant density ISM, Weaver et al. (1977) show 
that, even during the adiabatic evolution of the system, when the 
fluid in region (c) does not cool efficiently, the thickness of re- 
gion (c) is only about 0.14iis. Also the self-similar solutions of 
Ryu & Vishniac (1991) and Koo & McKee (1992), for a range of 
power law density profiles, Pb{R) oc R^" , and wind luminosity 
L{t) oc t^, the thickness of region (c) is only ~ 0.09-0.16i?s. We 
verify this in the context of our model parameters in Appendix [B] 
Further, explosions in an expanding universe generically have the 
swept matter concentrated in a thin shell (Bertschinger 1985; Voit 
1996). These results encourage the use of the thin shell approxima- 
tion even when the shocked halo/IGM gas cannot cool efficiently. 
We generally find that the shocked halo gas can indeed cool effi- 
ciently, but the shocked IGM gas need not do so, especially if it 
cannot mix efficiently with the metal enriched galactic wind mate- 
rial in region (b). 



and will only escape as its increasing internal pressure drives out 
the external medium. 

Our subsequent treatment of the dynamics of the galactic out- 
flow assumes a thin shell approximation, analogous to the treat- 
ment of interstellar bubbles driven by stellar winds (cf. Castor et 
al. 1975; Weaver et al. 1977; Ostriker & McKee 1988; Tegmark et 
al. 1993). In this picture (see panel (a) in Fig.[TJ the wind blown 
bubble, at some stage, has an onion-like structure with 4 concentric 
zones: (a) an innermost region consisting of the galactic wind blow- 
ing out (called "free wind") (b) a hot bubble of shocked wind gas 
(the galactic wind entering the surrounding halo/IGM gets shocked 
at an inner shock at radius say _Ri) (c) a thin dense shell of shocked 
IGM/halo gas separated from the shocked galactic wind by a con- 
tact discontinuity at Rc and (d) undisturbed halo/IGM gas outside 
an outer shock at radius Rs . In what follows we shall describe an 
outflow with the above structure as a "pressure driven outflow". 

While the above picture is valid for most of the time, in the ini- 
tial stages, the hot gas could be a filled bubble of radius Rs, with- 
out the innermost shock at R\ . In this case, the SNe in the galaxy 
would just feed this bubble directly with its hot ejecta. But sub- 
sequently, as this bubble expands, its pressure decreases, and then 



3.2 Modeling the outflow dynamics 

The evolution of the spherically symmetric outflows is governed by 
the following equations in the pressure driven case as in Fig. la (cf. 
Weaver et al. 1977; Ostriker & McKee 1988; Tegmark et al. 1993), 

A'Rs ^ ^ttRI {Pb - Po) _ rhs{Rs){Rs - vo{Rs)) 
dt^ ms{Rs) ms{Rs) 

GMjRs) 

' 

^iRs) = eATxRlpB{Rs)[Rs - vo{Rs)). (6) 

Here, the subscript s represents the thin shell variables while a 
subscript b represents the bubble variables. In particular, Rs is the 
physical radius of the shell and Pb the pressure inside the bubble. 
Further, ms{Rs) is the baryonic mass in the thin shell, M{Rs) 
the total gravitational mass within 7?^ and pb{Rs) is the baryonic 
density of the ambient medium into which the supper bubble is ex- 
panding. 

The first term in the right hand side of Eq. |5j represents the 
outward force due to the difference between the hot bubble pressure 
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{Pb) and the outside medium pressure (Po). The second term takes 
into account both the drag on the outflow due to mass swept up 
from the surrounding medium and the addition of momentum of 
this medium to the shell. The last term is deceleration of the shell 
due to gravity. Further, vq is the velocity of the fluid outside Rs- 
As the wind propagates it sweeps up mass from the surrounding 
medium. We assume that a fraction e of this swept up mass is added 
to the shell and rest of it is incorporated into the hot bubble. The 
resulting evolution of the shell mass rris is described by Eq. |6|. 

The pressure Pb and the thermal energy of the bubble, Eb, are 
related by, 



Eb 



(7) 



with the assumption that the adiabatic index of the gas 7 = 5/3 
The evolution of the energy in the bubble is in turn given by 

dEb 



dt 



= L{t) - A(t, T, Z) - 4n{RiRs - RiRi)Pb. 



(8) 



Energy is fed into the bubble from the supernovae explosions in the 
host halo and is lost due to cooling and the PdV work done during 
the expansion of the bubble. In Eq. l|7j, we have taken into account 
of the fact that the wind may enter the hot bubble through an inner 
shock at Ri . This leads to a smaller hot bubble volume compared 
to the models that do not consider the inner shock. Further since it 
takes a finite time, (Ri/v^), for the free wind to propagate to Ri, 
the effective energy input to the bubble in a time dt is, L(t)dt — 
Lo{te)dte, where te — t — Ri/v^ is the time at which the wind 
leaves the galaxy and Lo{te) is its mechanical luminosity generated 
from the SNe explosions in the galaxy. Here is the asymptotic 
velocity of the wind from the galaxy before it encounters the shock 
at Ri . As noted earlier, such a thermally driven wind can arise when 
the bubble pressure drops sufficiently below the ram pressure of 
the gas outflowing from the galaxy. The luminosity output in this 
situation can also be written as Lo(te) ~ Mm{te)v'l^/2, where 
is the rate at which mass is ejected by the galaxy. The evolution 
equation for Ri is obtained simply from the jump condition across 
the inner shock, assuming it to be strong: 



Pb = 



3 M„(t 



4 47ri?' 



Vw — -R] 



(9) 



Earlier semi-analytical studies of SNe driven outflows from galax- 
ies, do not take into account the evolution of the inner shock 
at Ri. Indeed for decelerating self-similar flows, Ri grows with 
time slower than Rs (cf. Weaver et al. 1977), and so its evolution 
changes Pb negligibly, and hence the evolution of Rs- However, 
we find that the outflows accelerate in the early stages of evolution, 
basically because of an increasing L{t) and a steeply decreasing 
Pb(R)- In this case, we will see below from naive scaling argu- 
ments that Ri can approach Rs- In addition, as discussed above 
when there is efficient cooling of the bubble Ri approaches Rs and 
the 'pressure driven flow' transits to a 'momentum driven flow'. 
Therefore, it is important to take account of the dynamics of Ri to 
determine the correct evolution of the outflow. 

The cooling of the bubble (represented by A{t, T, Z) in 
Eq. |8|) is due to Compton drag against the CMBR, bremsstrahlung 
and recombination line cooling. The cooling rate depends on the 
hot bubble density, temperature and metallicity (Z). The metallic- 
ity of the outflowing gas and interstellar medium of the galaxy are 
computed in a self-consistent way assuming instantaneous uniform 
mixing (see Appendix |A|. As expected the metallicity grows with 
time. However, to keep our computation simple, we estimate bubble 



cooling rate as a function of time for an assumed constant metal- 
licity (say Z — 0.01^0). We note that the maximum metallicity 
achieved in a typical bubble is of this order. Usually the bubble ma- 
terial reaches the maximum metallicity when it is outside the halo 
where the adiabatic cooling dominates over the radiative cooling. 
Thus assuming this constant metallicity slightly over estimates the 
cooling in the early stages of the wind and have no effect in the lat- 
ter stages. We adopt here the cooling rates A, given by Sutherland & 
Dopita (1993). The temperature and bubble density (required also 
for estimating A) are computed as follows. The mass of the hot 
bubble evolves as. 



M,„ dt 



+ / {l-e)4nR^spB{Rs){Rs 



dt. 



(10) 



Here, the mass outflow rate from the galaxy is A/„ which we as- 
sume to be proportional to the SFR ; i.e = tjMsf, where r] is 
the mass loading factor. Most outflow models that are discussed in 
the literature use 7; > 1 (see for example, Furlanetto & Loeb 2003, 
Oppenheimer & Dave 2006; Bertone, De Lucia & Thomas 2007). 
The value of rj is found to be greater than 1 for local star forming 
dwarf galaxies ^Martin I999[ >, and in the range 0.1 < rj < 0.7 for 
ultra luminous infrared galaxies at 0.04 < z < 0.27 l |Rupke et al.| 
2002). Martin (2005) finds a median values of r; of 0.19 and 0.09 
for high and low redshift ultra luminous infrared galaxies respec- 
tively. Even though there are predictions of the mass dependence of 
rj (see Murray, Quataert & Thompson 2005) available observations 
do not strongly support such a dependence l |Rupke et al. 2002^ In 
our model calculations we assume -q to be independent of mass and 
use T) — 0.3 in most of the models below. 

The second term in Eq. l |IO[ > takes into account the mass load- 
ing of the hot bubble due to processes like, evaporation of the 
shell, and we have assumed simply that the fraction (1 — e) of the 
halo/IGM gas that has not being swept by the shell is added to the 
hot bubble (cf. Furlanetto & Loeb 2003). The hot bubble is taken to 
be of nearly uniform density given by Pb = 3Mb/i47v{Rl-Rl)y, 
the self-similar solution derived by Weaver et al. (1977) for ex- 
ample, supports such a view, with only the region near the contact 
discontinuity between regions (b) and (c), being at a higher density. 
Knowing the bubble density, we can also derive its temperature Tb, 
from the total thermal energy, using 3k{Mb /'mp)Tb — Eb where 
m.p is the proton mass. Here we have assumed the gas is mostly 
ionized hydrogen. 

Given the temperature, density and metallicity of the bubble 
gas, we can calculate its cooling rate. We show below that for cer- 
tain range of parameters, the bubble can cool efficiently enough 
that Ri becomes very close to Rs- The outflow then transits to a 
'momentum driven flow'. The wind then directly deposits momen- 
tum at the shell and the evolution of Rs is governed by a modified 
equation (Ostriker & McKee 1988; Bertone, Stoehr & White 2005) 

d^Rs _ M^{t)v^{l - Rs/v^) _ 4-kRIPo 
dt'^ ms{Rs) irisiRs) 

ms{Rs)iRs - vo{Rs)) GM{Rs) 



ms(Rs) 

where, the shell mass m. now evolves as. 



Ri 



dms 
~df 



{Rs) = + 4-kRIpb{Rs){Rs - vo{Rs)). 



(11) 



(12) 
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We always begin the evolution of outflows using Eq. |5j and Eq. |6j 
and switch over to Eqs. ijTTJ and ljT2j, if and when R\ becomes very 
close to Rs . 

We also have to specify the distribution of the dark matter and 
baryons in the halo through which the outflow propagates. The dark 
matter distribution within the virial radius of the halo is assumed 
to be a NFW density profile ( [Navarro, Frenk & White \991) and 
smooth outside the halo with mean cosmological density. The bary- 
onic density ps is estimated in the following manner. Within the 
virial radius, a fraction fh ~ 0.1 of the total baryonic mass is taken 
to be still in the halo in hydrostatic equilibrium with the dark mat- 
ter potential at the virial temperature Tvir. This can represent for 
example, the gas which is being accreted into the halo in the hot- 
mode of accretion (cf. Keres et al. 2005; Dekel & Bimboim 2006) 
and which has not yet cooled and fallen into the galaxy. We note in 
passing that earlier work on outflows have adopted values ranging 
from fh = 1 (Madau, Ferrara & Rees 2001) to /h = (Furlan- 
etto & Loeb 2003). Kobayashi, Springel & White (2006) have used 
fh = 0.1 with an NFW profile in their isolated disk models. A gas 
in hydrostatic equilibrium in an NFW dark halo with virial radius 
7?vir and concentration c, has a density profile that is well fitted by 
a beta model (Makino, Sasaki & Suto 1998): 



Pb{R) = 



[1 + {R/R. 



(13) 



where, for a gas at the virial temperature, the core radius is Rc — 
(0.22/c)7?vir and pc is the central gas density, determined from 
normalizing the gas mass in this profile to the total mass. We have 
taken a typical concentration parameter of 4.8 for all the halos 
(Madau, Ferrara & Rees 2001). Hence ps is equal to the left over 
gas density in the halo for Rs < Rvh and outside Tivir, it is simply 
the background IGM gas density. Note that at i?vir the gas density 
is still larger than the background IGM density. To avoid any un- 
physical jump we assume an exponential decay of density of width 
0.2_Rvir (Madau, Ferrara & Rees 2001). In what follows we will 
examine the sensitivity of our results to fh ■ 

The outside pressure Pq is therefore fixed to be the pressure of 
the halo gas at the virial temperature (Tvir) within the virial radius. 
Outside this radius Po is calculated assuming that the gas in the 
IGM is at 10* K. The latter assumption is justified as the ioniza- 
tion front from the galaxy moves faster than the outflow and hence 
the wind always passes through an ionized medium which has also 
been photoheated to a temperature of 10* K. 

Further, in Eq. |6|, vo is the velocity of the surrounding 
medium. We adopt the following form for vo (Furlanetto and Loeb 
2001), 



vo{R) = 



0.0 for R < i?vi 
a /' R 



(14) 



4 for i?vir < -R < 4i?vir 



= 1^H{R - 4i?vir) for 4J?vir <R< 12J?vir 

= HR for 12i?vir > R. 

Thus the medium is assumed to be at rest within the halo. 

The mechanical luminosity, Lo{t), fed into the wind com- 
ing from the SNe produced by a continuous star formation (as in 
Eq. ^) is given by, 

fit, 



Lo{t) = 10' 



exp 



X u ft 
t{z) — t{Zc 



-M 



t{z) - t{Zc) 



erg s 



(15) 



Here, we have assume that each SNe produces 10^* ergs of energy 
and a fraction ew of this energy goes to power the wind. For most 
of our work we take — 0.1. Note that, Mori et al. (2002) from 
numerical simulation find an efficiency of 20-30% of converting 
supernova energy into kinetic energy of the outflowing gas. This 
is factor 2-3 higher than the maximum we use in this work. 
Further, is the number of SNe per unit solar mass of stars formed. 
For a Salpeter IMF with mass range 1 — 100 Mq, one SNe occurs 
every jy"^ ~ SOA/q of stars formed. However, if we assume 
Salpeter IMF with mass range 0.1-100 A/q then i/~* ^ ISOAfg. 

For the major part of the present work we adopt the contin- 
uous star formation model above, with the parameters constrained 
by fitting the observed high redshift luminosity functions (see Paper 
I). However, much of the earlier semi-analytic models of outflows 
assume that the star formation in a halo occurs in a single instan- 
taneous burst ( [Scannapieco 2005] [Madau, Ferrara, & Rees 2001| l. 
To compare our results with previous work we will also show the 
results obtained with such a model, where we assume that a frac- 
tion esF of total baryonic mass goes into stars instantaneously. In 
such a model, the last SNe explosion will occur at a time ~ 33 
Myr (assuming that lower mass limit of a star to explore as SNe is 
~ SMq) after the burst of star formation. This gives rise to 

Lo{t) = 10 X Ew esFT^Af^- ^ergs (16) 

"m lOB 

where toB is the maximum life time of an OB star which is as- 
sumed to be 33 Myr. 

Finally to solve the above listed equations one has to specify 
the initial conditions. 

• The initial radius of the hot bubble is taken to be Ri — 
J?vir/15. This is much smaller than halo size but typically larger 
than the radius of a disk galaxy forming in the halo (Efstathiou 
2000). 

• We wait until the bubble pressure determined from Lo{t) is 
equal to the outside thermal pressure of the halo gas. This time 
is taken to be the initial time ti for the integration of the outflow 
equations. Typically ti <^ tdyn and so most of the star formation 
and energy input into the hot bubble will take place at later times. 

• The initial mass loaded into the bubble is taken to be rj times 
the mass of formed stars upto ti, motivated by the proportionality 
of the wind mass loss rate and the SFR noted above. 

• The initial mass of the shell is taken to be equal to that of the 
super bubble at ti. However, we find that the subsequent swept up 
mass by the shell form the ambient medium greatly exceeds this 
initial mass (even after the first few time-steps), and so the evolu- 
tion of the outflow is almost independent of this initial assumption. 

• Although the shock at Ri arises perhaps after an initial period 
of filled bubble evolution, we shall solve the equations above as if 
it was always present. This only changes the initial period of the 
evolution. 

• We start the evolution with a initial zero velocity for the shell. 
Again, we will show that starting from a supersonic shell velocity 
makes very little difference to the final evolution (mainly because 
the bubble pressure at ti is not large enough to maintain this veloc- 
ity). In passing we note that even for subsonic expansion of buoyant 
bubbles into uniform fluid, one gets an evolution equation very sim- 
ilar to the equations we follow here (see Batchelor, 1997, p. 479), 
including an effective mass loading due to the displaced mass. We 
also set Ri — initially. 

We follow the evolution equations upto when the peculiar ve- 
locity of the shock reaches to the local sound speed and afterward 
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Figure 2. Properties of an outflow originating from a galaxy of dark matter 
halo of mass IO^Mq that has collapsed at = 6. Panel (a) shows the time 
evolution of the shock velocity (Vs), while panel (c) shows the outer shock 
location (Rg ) and the inner shock location ) as a function of time. The 
velocity Vs and the ratio Ri / Rs as a function of the outer shock location 
are shown in panel (b) and (d) respectively. The vertical lines in panel (a) 
and (c) mark the dynamical time-scale for this halo. The vertical lines in 
panels (b) and (d) and the horizontal line in panel (c) mark the virial radius. 

allow it expand with the Hubble flow. Note that if the velocity of 
an outflow becomes zero within the virial radius of the halo and 
the elapsed time is less than tdyn, we start the outflow again from 
Ri. The reason is that the Lo{t) increases upto tdyn in our model 
and hence there is stifl possibihty of creating a new supper bubble 
which can escape the halo potential. 

Having drawn the basic framework of our calculations, in the 
following section we discuss various generic properties of the out- 
flow solutions in our models. We also test our model equations, 
and its starting conditions by comparing to the known self-similar 
evolution when both Lo ) and pb{R) have a scale-free power law 
form. We show that the outflow dynamics is robust to fairly large 
variations of Ri, Vi and ti around the above fiducial values. 



4 STRUCTURAL PROPERTIES OF THE OUTFLOW 

In this section, we illustrate typical structural properties of galactic 
outflows, by focusing on the case of a 10^ M© halo collapsing at 
Zc = 6. We will find that halos having masses around this value, 
dominate in determining the global consequences of outflows. For 
this halo the virial radius is i?vir ~ 4.7 kpc. The evolution of var- 
ious physical quantities are shown in Fig. |2] and Fig. |3] We have 
taken /, = 0.5, k = 1, fh = 0.1, e = 0.9, = 0.1, r? = 0.3, 

= 50 Mq and the initial conditions are set as described above, 
with Ri = i?vir/15, Vi = 0. We call it as our fiducial model. 

Panel (a) of Fig.|2]shows the time evolution of the shell veloc- 
ity (lis). The outflow gets accelerated initially, due to the increasing 
energy input from the galaxy and decreasing halo density profile. 
This phase lasts typically for a timescale ~ t^yn when Lo{t) in- 
creases, after which the outflow starts to decelerate. Further, the 
acceleration period is mostly when the shock is within the virial ra- 
dius of the halo or very close to the virial radius. To make this more 
clear we have plotted the velocity of the outflow as a function of 
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Figure 3. The time evolution of hot bubble temperature (solid line) and 
shell temperature (dashed line) are shown in panel (a) for the outflow origi- 
nating from 10^ A/0 halo that has collapsed aX Zc = 6 

. Panel (b) shows the metallicity evolution of the bubble for the same out- 
flow. The vertical dotted lines in both the panels represent the dynamical 
time-scale for the halo. 

location of the outer shock in panel (b) of Fig.|2] The virial radius 
is shown here with the dotted vertical line. At the virial radius the 
density of outside medium changes sharply, which results in steep- 
ening in the slope of the velocity profile. The outflow subsequently 
decelerates till the peculiar velocity decreases to the sound speed of 
the IGM, at which stage we assume it freezes into the Hubble flow. 
In the above example this happens at 3 x 10^ years when the radius 
of the outflow is ~ 290 kpc. By z = 1 the outflowing material has 
spread to a proper distance of ~ 370 kpc from the galaxy. 

In panel (c) of Fig.|2]we show the time evolution of Ri (dashed 
line) and Rs (solid line). In this model Ri becomes ~ O.GRs over 
a period of few dynamical time scale, as seen from panel (d). When 
the injection of mechanical luminosity stops Ri becomes zero as 
there is no ram pressure to support the bubble pressure. 

The evolution of the bubble temperature, Tt, (solid line in 
panel (a) of Fig. |3]l is governed by competing effects. Adiabatic 
expansion and radiative cooling lead to a decrease of T^, while en- 
ergy input from the galaxy can result in its increase. The radiative 
cooling depends on the metallicity of the gas. In panel (b) of Fig.|3] 
we show our computed bubble metallicity as a function of time. 
The details of our metallicity calculations are given in Appendix[A| 
In most of our models we use the Carbon as the metallicity in- 
dicator with an yield of 0.1 Mq per solar mass of SNe. Initially 
the metallicity increases with time as most of the bubble material 
comes from the galaxy. The maximum metallicity of the bubble 
gas is ~ 0.05 Zq. In the later stages the metallicity begins to de- 
crease when the star formation in the halo stops while more and 
more primordial gas gets added to the bubble. In the calculations 
presented above, the radiative cooling of the bubble gas is calcu- 
lated assuming the average metallicity of the gas to be 0.01 Zq. 
In Fig. [3] Tb starts from an initial value ~ 5 x lO*" K (determined 
by the mass loading factor 77), decreases for t ~ 10^ yrs, before 
increasing again due to the growing SFR of the galaxy till t ~ idyn 
(Eq.[2|. After t ~ tdyn the SFR starts to decrease and the effect of 
adiabatic expansion dominates over the energy input to the bubble, 
leading to decrease of Tt. We notice that the results do not change 
even when we use the average metallicity to be 0.1 Zq. This is 
mainly because at T;, > lO'' K the adiabatic cooling is much faster 
than the radiative cooling. 

The dashed curve in panel (a) of Fig.[3]is the post shock tem- 
perature of the shell assuming it to be an adiabatic shock. As this 
temperature depends on the velocity of the shell it follows the time 
evolution of the Vs shown in panel (a) of Fig. |2] In the case of thin 
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Figure 4. Top and bottom panels show the cooUng time for the shell origi- 
nating from halos of mass 10^ Mq and 10^^ Mq respectively. The solid 
and dashed curves are for primordial abundance and Z = 10~^ Zq re- 
spectively. The dotted lines in both the panels correspond to cooling time 
equals to the expansion time. The collapse redshifts, Zc, are indicated in 
each panel. 

shell approximation we expect the swept up shell material to cool 
very efficiently. We check this in the following subsection. 

We note that the pressure of the bubble Pj, is always a mono- 
tonically decreasing function of time. This arises due to the fact that 
the density inside the bubble always decreases and this decrease is 
faster than the increase of temperature at any time in the evolution 
of the outflow. 



4.1 Cooling of the shell gas 

It is also of interest to examine whether the swept up gas in the shell 
can cool efficiently. In Fig.|4]we show the cooling time (tcooi) of the 
shell gas as a function of the expansion time (toxp) of the outflow, 
assuming the shell density is simply 4 times the pre-shock density, 
(as would obtain for a strong adiabatic shock of a 7 — 5/3 gas) and 
its temperature is the post-shock temperature. The cooling rate has 
been taken from Sutherland & Dopita (1993) assuming primordial 
abundance (solid lines) or Z = 10~^ Zq (dashed lines). The latter 
case is to illustrate the effect of having some metals in the halo 
gas, perhaps due to enrichment from earlier generation of outflows, 
or from partial mixing with the metal rich bubble. We have also 
plotted as a dotted line the relation tcooi = tcxp. The top panel 
is for a 10® M0 halo that has been discussed in detail above. The 
bottom panel is for a IO^^A/q halo collapsed at = 4 with all 
other parameters same as that of the fiducial model used for the 
illustration above. 

One sees that in the absence of radiative heating by the UV 
flux from galaxy the shell can indeed cool efficiently with tcooi < 
texp while the shock traverses much of the halo. The epoch when 
the shock leaves the halo is the time when tcooi rises abruptly. Just 
around this epoch or just before, one gets fcooi ~ tcxp. One can also 
see that the cooling becomes inefficient in the IGM (tcooi > toxp). 
Note that the density in the shell is not expected to be uniform 
while the outflow traverses the declining density distribution in the 
halo. Indeed self-similar solutions for the adiabatic evolution of 
the shell in a medium with declining density profiles, by Ryu and 
Vishniac (1991) show that the density (temperature) in the shell 



rises (falls) steeply as one goes from the shock at Rs to the con- 
tact discontinuity Rc- For example, for the case L{t) = constant, 
Pb{R) oc (see Fig. le in Ryu and Vishniac 1991), the den- 

sity (temperature) is a factor of 2 larger (smaller) than the value at 
the shock, for the inner half of the shell. We have also obtained the 
shell structure assuming an adiabatic strong shock for our model 
parameters in a self-similar outflow (see appendix |BJ. When the 
outflow is traversing the halo, we also get a similar steep density 
enhancement in the shell compared to its density at the shock. As- 
suming A oc T~^''^ in the relevant temperature range, the den- 
sity and pressure profile together make the icooi smaller by a factor 
^ 5.7 for the inner half of the shell. So our estimates of the shell 
cooling times, given in Fig.|4] are more upper limits while the out- 
flow traverses the halo. In addition, since the shell/hot bubble in- 
terface is likely to be unstable, leading to mixing (see below), the 
metallicity of the shell gas could also be much higher than we have 
assumed. Overall it appears reasonable to assume that the shell can 
cool efficiently while it is traversing the halo, but not in the IGM. 
Later when the outflow has gathered up mass from the IGM and 
slowed down sufficiently, tcooi can again become smaller than tcxp. 
In case of efficient cooling, the thin shell density will be enhanced 
to ps ~ pbM^ , where M is the Mach number of the shock. These 
results will be useful in our discussions below. 

4.2 Transition from pressure to momentum driven flow 

We note that the basic trend seen in the illustrative example dis- 
cussed till now is the case with most of our models. That is through- 
out the evolution, the flow is pressure driven. However, in few 
cases, R\ can become close to Rs and the outflow can transit to 
a "momentum driven case". This arises when the radiative cooling 
of the hot bubble gas is efficient. When the bubble is well within the 
halo its density and metallicity depends mostly on the wind mate- 
rial. As we have already mentioned, the recombination line cooling 
depends on metallicity, temperature and density of the hot bubble. 
Higher density and lower temperature that favors cooling can be 
achieved by increasing the mass loading of the bubble (i.e higher 
value of rf). 

To illustrate these points we have chosen three examples and 
show the time evolution of R\ and Ra in Fig. |5] For all the three 
examples, we have assumed /* = 0.25, k = 0.5, r) — 2.5, fh ~ 
0.02, e„ = 0.15, the metallicity of the bubble gas Z = 0.03 Zq 
and Zc = 4. The halo masses are indicated in the corresponding 
panel. The initial temperature of the gas is ~ 10'^ K in all the mod- 
els as 77 is same in all cases. However, bubble density is higher in the 
case of higher mass halo. In the first case, for M = 5 x 10^" Mq, 
the flow is completely pressure driven (panel (a) of Fig.|5]l. Here 
Ri is ~ 0.7iis. The time evolution of various physical quantities 
follow the example discussed before. In the second example, for 
M = 10" Mq, 7?i ~ 7?3 at t = 10* yr (panel (b) of Fig.jl}. Af- 
terward the flow transits to a momentum driven wind case. In our 
models the mass loading into the bubble is proportional to the halo 
mass. The main difference in this model compared to the model 
discussed in panel (a) is that the bubble density is a factor 4 higher 
in the second case. This makes the radiative cooling rate a factor 4 
times faster. This is sufficient to cool the bubble gas to T = 10* K 
in ~ 8 X 10^ yrs. As the shell-l-bubble gas is already been acceler- 
ated to high velocity and mechanical luminosity is still raising (as 
t < tdyn) the gas continuous to move outwards. When we con- 
sider M = 5 X 10" Mq (panel (c) of Fig.js) the initial gas den- 
sity is higher by a factor 7. At this density the gas cools to 10* K 
within 1.5 X 10^ yrs. This time is very short and gas does not have 
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Figure S. Time evolution of Rg (solid lines) and R\ (dashed lines) for three different halo masses. All the curves are drawn for /, = 0.25, k = 0.5, 77 = 
2.5, fh = 0.02, = 0.15, bubble metallicity Z = 0.03 Zq and Zc = 4 but different halo masses as indicated in the corresponding panels. In the case of 
Af = 5 X 10^'' Mq the flow is a pressure driven flow through out (panel (a)). However for M = 10^^ Mq the flow transits from a pressure driven flow to 
the momentum driven flow (panel (b)). In the case of Af = 5 X lO'^^ Mq the outflow is not set as the gas cools much faster. 



sufficient energy to initiate an outflow. So the bubble gas remain 
confined to the galaxy itself. 

In principle one can also get pressure driven flow transiting 
into a momentum driven flow in cases where there is high mass 
loading with higher metallicity. However, we find that for a wide 
range of model parameters the pressure driven flow is more generic 
than a flow which transits to momentum driven case (especially be- 
cause we take 77 = 0.3 for most models). In any case, as shown in 
the above example, our calculations switch to correct set of equa- 
tions as and when the outflow transits to a momentum driven case. 



4.3 Comparison with scale-free solution 

The thin shell equations allow for a power-law solution when Lq (i) 
as well as Pb{R) are scale-free, provided we ignore gravity, out- 
side pressure and the effect the inner shock. Suppose the halo 
gas density goes as pb oc R~°' and Lo{t) oc t'^, then simple 
scaling argument suggests that the outflow radius will scale as 
Rs{t) oc {Lot^/pBY'^ oc r, where n = (3 + /3)/(5 - a). 
Acceleration of the outflow obtains for n > lora + /9 > 2, 
while deceleration obtains in the opposite limit. For example, if 
Q = 2.8 and /3 = 1 (as one expects within the virial radius of the 
halo and at early times in our model), one expects Rs{t) oc t^^^'^, 
Vs{t) = Rs{t) oc andP6(f) = Eb/{2i:Rl) oc t-^^/". 

As a simple test of some aspects of our model, we first show in 
Fig. [6] the evolution of an outflow with I/o(i) oc t and pB oc R~^'^ 
and compare it with the expected scalings of Rs (t) , Vs (t) and Ft (t) 
derived above. Here we have also used the initial conditions as ex- 
plained earlier. We see that, after a short initial period, there is in- 
deed excellent agreement of the computed evolution of Rs , Vs and 
Pb, with the expected scaling laws derived above. 

We also examine in Fig. |7] the actual scaling behavior of the 
outflow radius, Rs{t), for a halo of 10^ Mq considered above, 
now incorporating gravity, external pressure and the effect of the 
inner shock. One sees that, inside Rvir, n is somewhat smaller 
than that expected for the pure self-similar evolution law, assum- 
ing pB oc R^^'^ and Lo{t) oc t. This is to be expected due to the 
influence of gravity and outside pressure, which both go to slow 
down the expansion of the outflow. At distances much larger than 
7?vir, the gas density decreases with the expansion of the universe 



approximately as pB oc and Lo(i) = 0. Even if we assume 
negligible radiative cooling of the hot bubble, Et would decrease 
due to the expansion of the universe. This case has been analyzed 
by Voit (1996), by a conformal transformation of the fluid equa- 
tions. Voit finds that in the scaled variables the structure of the 
shocked IGM gas would follow the usual Sedov self-similar so- 
lution. In terms of the actual time variable, Voit predicts that for a 
flat matter dominated universe, Rs{t) oc [1 - (t/uy^^^f^^t^^^ . 
Such an evolution is indeed close to that obtained for the actual 
outflow, as can be seen from Fig.|7] 

We now examine the sensitivity of the outflow evolution to the 
initial conditions that we have adopted. 



4.4 Sensitivity to tlie adopted initial conditions 

Let us begin with the sensitivity of the outflow solution to the 
adopted initial radius of the supper bubble. In panel (a) of Fig. [8] 
we show the shock velocity Vs = iis as a function of the shock 
radius Rs(t), for a 10^ Mq halo collapsing at Zc — 6 discussed 
above. The solid line is for Ri = _Rvir/15 while the dashed lines 
is for Ri = _Rvir/5 keeping all other parameters as in the fiducial 
model. It is evident from the figure that a modest change in Ri has 
very little effect in the evolution of the outflow. 

Panel (b) of Fig. [8] looks at the effect of changing the initial 
velocity Vi from Vi — (solid line) to Vi — 3uc (dashed line) with 
other parameters as in the fiducial model. Starting the evolution 
of the outflow with such a larger Vi has very little effect on the 
time evolution of the outflow. This is because the mass loading in 
the shell very rapidly decreases the velocity to that which can be 
consistently maintained with the existing bubble pressure. 

Finally, panel (c) of Fig.[8]shows the effect increasing ti by a 
factor 5. After an initial period the shell velocity again latches on 
the fiducial model solution. And the radius to which the wind prop- 
agates is almost independent of ti for a reasonable range around 
our fiducial value. Therefore we conclude that the final results of 
our models are almost insensitive to the initial conditions, for rea- 
sonable variations around our fiducial values. 

In the following section, we study how the model parameters 
influence the the outflow properties, in particular the extent of the 
outflow (Ra). 
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Figure 8. Effect of using different initial conditions on tlie outflow evolution. Panel (a) shows the effect of changing the initial radius from iJvir/15 to i?,vir/5. 
Panel (b) shows the effect of changing the initial velocity from to 3dc. Panel (c) shows the effect of changing initial starting time. In all the panels the results 
for our fiducial model are shown by the solid lines and the evolution of the modified models are shown by dashed lines. 
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Figure 6. The asymptotic self-similar behavior of properties of the outflow. 
The position of outer shock (Rs), velocity (Vs) and hot bubble pressure 
(Pfj) of the outflow coming out from a halo of mass lO'' Mq are shown in 
bottom, middle and top panel respectively. The input energy is assume to be 
Lo{t) oc t and the outside gas density is oc R~^'^. The corresponding 
expected self-similar power-law solutions are shown by the dashed line. 
The power-law relationships expected from the scale free solutions are also 
given in the corresponding panels. 



5 DEPENDENCE OF OUTFLOW PROPERTIES ON 
MODEL PARAMETERS 

The free parameters in our model are, k, IMF, e, rj, fh, and £„, 
in addition to the background cosmological parameters. As pointed 
out before, for most part of this paper we will use cosmological 
parameters from the WMAP 3rd year data and /, — 0.5 and k = 1 
that fits the observed galaxy luminosity functions at 3 < z < 6. We 
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lO"^ 10^ lO'^ 

Time (yr) 

Figure 7. The actual solution for the time evolution of Rs in our model. 
The expected self-similar behavior are shown for two different regimes. 
Within the halo, Rs{t) oc t^l^-'^ (short dashed line) and far away from 
the halo, Rs(t) oc [1 - itltiY^I'^^fl-n^l'-^ (long dashed line) with ti = 
5 X 10'^ yrs. 

are left with IMF and 4 parameters of the model that are associated 
with wind dynamics. In this section we explore the dependence of 
our results on the choice of these parameters. In addition we also 
explore models that have star formation in a burst mode. 

5.1 The halo mass fraction 

We first consider how the outflow properties change when we 
change the halo fraction, jh- Fig. |9] shows the velocity evolu- 
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Figure 9. Effect of halo mass fraction on the outflow. The top panels show 
the evolution of Vs as a function of time (left panel) and Rs (right panel) of 
the outflow coming out from a halo of mass 10^ Mq . The solid lines are for 
fh = 0.1 and the dashed lines are for = 0.3. In the bottom panels we 
show the corresponding curves for the outflow origination from 10^^ Mq 
halo. 



200 


1 ^ a J 












100 













1 10 100 
{ kpc ) 




10 100 



lO'^M-, ,z -4 





Time ( yr ) 



R, ( kpc ) 



Figure 10. Effect of r) on the outflows. Panel (a) and (b) show respectively 
the shell velocity, Vs and hot bubble temperature, Tj, as a function of Rs 
for the outflow coming out from a 10^ Mq halo. Panel (c) shows the time 
evolution of Vs while panel (d) shows the Vs as a function of Rs for the 
outflow originating from a 10^^ A/q halo. In all the panels solid lines are 
for 77 = 0.3 and the dash lines are for r; = 1.0. 



tion characteristics for two halos of masses 10 (top panels) and 
10^^ Mq (bottom panels), for both fu =0.1 (solid hne) and 
fh ~ 0.3 (dashed lines). Higher value of fh increases the pres- 
sure in the halo (Po)- Thus our initial condition that bubble pressure 
equals to the halo pressure is achieved in a latter time compared to 
the model with smaller fh- The mass loading of the shell when it 
propagates inside the halo is larger in the case of high fh ■ For the 
IO^Mq halo, this leads to a lower peak velocity in the case of high 
fh (see upper panels of Fig.|9]l. It is however interesting to see the 
final velocity and radius of the outflow remains the same in both 
cases. This is mainly because when the shell travels in the IGM 
most of its mass comes from the swept up IGM material and the 
final velocity and radius are insensitive to fh- The above conclu- 
sion holds provided that the initial decrease of Vs does not lead to 
confinement of the outflow within iivir- Such a confinement is in 
fact obtained for the higher mass halo shown in the bottom panels 
of Fig|9] For this halo of mass 10^^ Mq, fh = 0.1 leads to an out- 
flow which escapes the halo while increasing to an = 0.3 leads 
to confinement of the wind. Hence the value of could crucially 
impact upon whether one has an outflow or not, especially for high 
mass halos. 



5.2 Mass loading from the galaxy 

Here we discuss the dependence of the outflow evolution on the 
mass loading factor, 77. The top panels of Fig. |10| show the results 
for M = 10® Mq for ri = 0.3 (continuous curves) and 77 = 1.0 
(dashed curves). We do not notice much difference in the evolu- 
tion of Rs however the temperature of the bubble in the case of 
77 = 1 is lower than that in the case of 77 = 0.3. Higher values of 77 
reduces the initial temperature and increases the density of the bub- 
ble when it is inside the virial radius. The bottom panels in Fig.|10| 



shows the results for M = 10^^ Mq- When 77 = 0.3 the model 
produce outflow that escapes the galaxy. However, for 77 = 1.0 the 
flow is confined. This is mainly because in addition to the initial 
temperature being low the bubble gas cools much more efficiently 
in the case of 77 = 1.0, due to its increased density. The other quan- 
tity that changes appreciably to the change of 77 is the asymptotic 
metallicity limit. This is higher in the case of high 77 (see Eq. (|A3|). 
Like in the case of fh, the influence of 77 on the nature of the flow 
seems to be strong when one considers high mass halos. 



5.3 Energy input efficiency 

Now consider the effect of varying the efficiency of energy input 
into the wind which is characterized by the factor (e^i/f,)- This 
can arise for example in changes to the adopted IMF, since such a 
change will lead to different v- Al 100 AIq Salpeter IMF gives 
I' — 1/50, while changing the lower mass cut-off to 0.5 AIq and 
0.1 A/q leads to = 1/70 and v — 1/130 respectively. In most 
of our models discussed till now we have assumed to be 0.1. 
A smaller value of this efficiency factor will have an effect similar 
to reducing the lower mass cut-off in the IMF. In the top panel of 
we have shown the results from a halo of mass o 
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Fig. 

with the above three IMFs. Two major effects are noticeable. The 
maximum velocity attained by the outflow and the final outflow ra- 
dius are both higher for higher 1/ (i.e IMF with higher values of low 
mass cut off), as expected because of the increase in input energy. 
Note that the same effect can be obtained by lowering the /, (or em) 
by a factor 2.6. In the case of 10^^ Mq (bottom pane of Fig. [Tl} 
the material is confined to the halo for f = 1/130. Remember 
this corresponds to standard Salpeter mass function. In Paper I, we 
have argued that the redshift evolution of UV luminosity function 
of galaxies can be explained by slowly evolving low mass cut of in 
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Figure 11. The effect of IMF on outflows. The top panels show (panel 
(a)) and Va (panel (b)) of an outflow coming out from a halo of mass 
10^ Mq collapsed at 2c = 6. The solid, dashed and dotted-dashed lines are 
for 1 - 100 Mq, 0.5 - 100 Mq and 0.1 - 100 Mq Salpeter IMF respec- 
tively. The bottom panels show the con'esponding results for the outflow 
coming out from a 10^^ Mq halo collapsed at 2c = 4. 



the IMF over the redshift range 3 to 6. Such an evolution will have 
important effects on the influence of winds in the global properties. 



5.4 Burst versus continuous made of star formation 

We now consider the effect of burst mode of star formation on the 
wind dynamics. In Fig. [12] we show the evolution of an outflow 
emerging from a halo of mass 10^ Mq collapsed at Zc — 6. Panel 
(a) of Fig. 1 12| shows the radius, Ks of the outflow as a function of 
time for esF ~ 0.5 (short-dashed line), 0.25 (long-dashed line) 
and 0.1 (dotted-dashed line). The other parameters are same as our 
fiducial model. For comparison we have shown as a solid line, the 
evolution of the same outflow in a continuous star formation model 
with /* = 0.5. It is evident from the figure that the maximum 
velocity achieved by the outflow is much higher in the case of burst 
mode of star formation compared to the continuous one. However 
the final radius of the outflow is smaller in the burst model. Further, 
as one expects, both the maximum velocity and the final radius of 
the outflow decrease with lower values of esp- Also note that the 
duration of initial acceleration phase is shorter in case of the burst 
model. The smaller size to which burst models drive outflows will 
have consequence on their volume filling factor (see below). 



5.5 Results for different halos 

We now consider outflows originating in halos of different mass 
scales. We assume the parameters of SFR and wind models iden- 
tical to the fiducial model discussed in section |4] For each mass 
we take the collapse redshift to be approximately the epoch when 
a 3(T fluctuation becomes unity. In Fig.[T3]we have shown the time 
evolution of the velocity (vs) and radius (Rs) of outflows arising 
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Figure 12. The evolution of an outflow originating from a 10^ Mq halo 
for the burst mode of star formation. The time evolution of Ks is shown in 
panel (a) while panel (b) shows Va as a function of Ks ■ The short-dashed, 
long-dashed and dotted-dashed lines are for esF = 0.5, 0.25 and 0.1 re- 
spectively. The other parameters are same as our fiducial model. For com- 
parison we have shown as a solid line, the evolution of the same outflow in 
a continuous star formation model with /* = 0.5. 



400 





400 




400 



200 



M = lo' M,, 

- 1 iiiiiii 1 iiiiiii TrniipT™ 


M^lO'^M., 
" z^ = 8 ^ 

: /\ ' 


7 M= IO'm. 

" Zj = 6 7 
A 

< 1 ^ ^ - 


- 1 iiiiiii] 1 1 1 llllfTi7liHl[ 1 Mill 
7 M = lo'°M. 
" z,, = 5 

/ A 
/ / 


."1 1 mill 1 1 iiiiTipriiiiiifTiira 

" M= lo" M , , r 
'\ /- 

/ 

^ y / - 


- 1 iiiiiii] 1 iiiinfMiiiy[ i m 

" z =3 V " 
/ 

/ / - 
/ / 



800 



400 



800 



Pi 



10' 10° 10' 10' 10° 10^ 

Time ( yr ) 

Figure 13. The time evolution of the radius (solid lines) and the velocity 
(dashed lines) of the shell arising from halos of different masses. The mass 
of the halo and its collapsed redshift are given in each panel The horizontal 
dotted lines in each panel represent the virial radius of these halos. 

from halos of different masses, collapsing at redshifts, Zc, as in- 
dicated in the figure. All other model parameters are taken to be 
same as the fiducial model. It is evident from the figure that the 
over-all qualitative features of outflow properties remain the same 
for all the relevant halo masses (for the assume set of parameters). 
The final radius of the wind and also the maximum outflow veloc- 
ity increases with the mass of the halo. This is a manifestation of 
the fact that SFR hence the energy available for the wind is higher 
for higher mass object when the bubble cooling is not that efficient. 
The halo gravitational force which increases with mass could have 
reversed this trend, but is generally sub-dominant compared to the 
pressure of the hot bubble driving the outflow, provided the outflow 
escapes the halo. We also notice in the previous sections, changing 
the parameters from the values used in the fiducial models may stop 
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the large scale outflows in the high mass halos. Basic requirement 
in this case is that the bubble cooling should be fast enough. This 
can be achieved in the high mass halos by increasing rj or fh or by 
decreasing e-u,- We will fold in these properties of different mass 
halos to calculate the global effects of outflows below. 



6 RAYLEIGH-TAYLOR INSTABILITY AND SHELL 
FRAGMENTATION 

We have seen in the last section that when the large scale outflows 
are possible, the shell of gas in the outer shock in our models gener- 
ically accelerates till just beyond the virial radius. This arises be- 
cause the input kinetic luminosity increases with time, till t = tdyn, 
while at the same time the confining halo density decreases with 
radius. As mentioned earlier, acceleration of the outflow obtains 
for n > lora + /3 > 2 (see section [43^ . In our case where 
Q ~ 2.8 within the virial radius, and /3 ~ 1 for t < tdyn, one 
gets n ~ 1.8 and so the shell accelerates while it traverses the 
halo. In the rest frame of the shell, this acceleration corresponds 
to a pseudo-gravitational force pointing from the dense shell to the 
low density hot bubble. Such an accelerating dense shell driven by a 
low density hot bubble is subject to the Rayleigh-Taylor (RT) insta- 
bility. (Physically, it corresponds to the analogous case of a heavy 
fluid lying 'on top' of a light fluid in a gravitational field). This can 
lead to fingers of the shell material penetrating the hot gas, while 
bubbles of the hot fluid rise into the shell, resulting in turbulent 
mixing of the two fluids. If the mixing scale becomes of order the 
shell thickness the RT instability will also lead to the fragmenta- 
tion of the shell. We now examine in a simple manner, some of the 
consequences which could result from this instability. 

In an expanding bubble, the evolution of the Lagrangian per- 
turbation ((/!>) in the shell with co-moving wavenumber kh is given 
by (Pizzolato & Soker 2006) 



a 

2- + 
a a 



2 J 







(17) 



where (cf. Chandrasekhar,1981; Padmanabhan, 2002) 

a p. 



(18) 



Ps — pb 

coth{khh/a) + pb ' 

Here, g is the gravitational acceleration of the halo potential, a is 
the expansion scale factor define as a{t) — Rs{t)/ Rs{to) and h 
is the thickness of the shell. We have also assumed that the hot 
bubble is "thick", that is khRs/a 2> 1. We take the initial time 
to as the time when we start calculating the growth of 0. We have 
also included above the effect of viscosity which stabilizes the RT 
instability for large fc^, in a heuristic manner akin to Piriz et. al. 
(2006). The effective kinematic viscosity Ve is obtained using (cf. 
Piriz et. al. 2006) 

_ 2(pb + Ps) 



pb + ps 

where pb and ps are the dynamical viscosity of the bubble ma- 
terial and shell material respectively while pb and ps the corre- 
sponding densities. For a fully ionized gas the Spitzer viscosity 
v = p/p^ 6.5 X 10^^(T/10^)^-^(ni/cm-^)-^ cm^ s-\ where 
T is the temperature and Ui the ion density of the gas. We have 
already seen the shell gas can cool faster than the expansion time 
while it is within the virial radius, to a temperature Ts ~ 10^ K. 
The density of the shell is taken to be ps — pB, where the 
Mach number Ms = (va/cs), with Cs the sound speed correspond- 
ing to a temperature Ts . So the above is largely determined by 
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Figure 14. Thickness of the shell (solid lines) and the mixing length (dashed 
line) for RT instability as a function of the outer shock location, for Rs < 
-Rvir- Left panel is for 10® Mq halo collapsed at Zc = 6 and right panel is 
for 10^^ Mq halo collapsed at Zc = 4. 
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Figure 15. The growth of perturbations at different physical length scale. 
We show the (f>{t) / (j){to) at t = 10^, 1.5 X lO*^, 2 X lO'^ and 2.5 X 10*^ 
yrs from bottom to top for 10® Mq (left panel) halo. For lO^^ Mq (right 
panel) halo wehave plotted(/)(t)/<^(to ) att = 4x 10'', 6x10'^, 8x10'' 
and 1.1 X 10* yrs from bottom to top. 



the bubble temperature and the shell density (since Tb 2> Ts and 

Ps > Pb)- 

In the absence of viscosity small scale perturbations always 
grow faster than larger scale ones. However, viscosity damps the 
growth of small scale perturbation. Expansion on the other hand 
damps the growth of perturbations on all scales. The initial per- 
turbation amplitude is unknown, but one does not expect the con- 
tact discontinuity to be smooth because of the small-scale inhomo- 
geneities in the wind fluid and the outside halo gas. When pertur- 
bations go nonlinear, such that the displacement of the shell-bubble 
interface (Rc surface) becomes comparable to the thickness h of 
the shell, one expects the shell to fragment. One can estimate h 
-iTxRlhps 



using rris 



This is shown in Fig. 

^9 
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for two repre- 
sentative halos, one of lO^M© (left panel) and for lO^M© (right 
panel). For the IO^Mq halo the thickness grows from h ~ lOpc to 
h ~ lOOpc, while for a lO^M© halo, ft ~ 100 - 200 pc. 

We can solve Eq. to find the evolution of 0(4) as a 
function of kh- In Fig. |15| we show the resulting growth factor 
G = (j){t) / (f){to) as a function of the proper wavelength of the per- 
turbation, Ap = {2n/kh)a, for various times t. We see that there 
can be significant growth of perturbations, by over an order of mag- 
nitude, due to the RT instability. For a IO^Mq halo, the most un- 
stable modes have Ap = Z/ ~ 10 pc, while for a lO^^Af0 have one 
gets a much larger // ~ 0.7 — 2 kpc. So the shell will keep losing 
mass due to the RT instability with the scale of fragments being 
comparable to // if all scales have similar fractional perturbations. 
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There could of course be a range of scales which become nonlinear, 
because we do not know the initial spectrum of inhomogeneities. 

Further, one expects the RT fingers/bubbles to have a random 
terminal velocity Vf ~ (Oef fif)^^^ (''f- Dimonte et al 2005 and 
references therein). Here gie// ~ |p|+J?a. Taking typical values for 
Qeff and leff, one gets Vf ^ 18 km s^^ {g^ff/10~'^ cm s~^)^^^ 
{If /WO pc)^/^ For a large enough Reynolds number Re = 
(vflf/ve), this contact layer between the two fluids could also be- 
come turbulent. In this case the growth of the dominant bubbles be- 
come self-similar, the two fluids mix in a mixing layer growing in 
size as Imix ~ 0.05(?e//i^ ~ T5 pc ((/e///10~** cm s~^) (t/(5 x 
10^ yr)) (Dimonte et al 2005). A gain if Imix becomes of order of 
the shell thickness the shell will be fragmented with the additional 
effect that the shell gas and the hot bubble gas within Imix also 
turbulently mix with each other. We show in Fig. [14] the evolution 
of the scale Imix with time. We see that Imix ~ h well within the 
virial radius, and so the shell is likely to fragment in this case due 
to the RT instability. 

The shell fragments resulting from the RT instability will be 
moving in the hot bubble and so will suffer significant mass loss 
due to evaporation. The rate of evaporation is given by (Cowie & 
McKee 1977) 

^ 5/2 

Me„ = 4.1 X 10"' ' 



\10^ K 



) (lO^)^-^'^" ^''^ 



where a is the radius of the fragment. This correspond to an evap 
oration time scale t^v = nic/M^v given by, 

-5/2 / 



t^v = 2.3 X 10° 



0.1cm- 



Tb 
106 K 



V 100 pc 



yr (20) 



where we have taken the cloud mass to be rric 



(47r/3)psa^ In 

fact if the bubble-shell interface becomes turbulent, such evapora- 
tion could be enhanced due to an enhanced contact surface. For 
10^ Mq halos, where ~ 10 pc, one expects the shell fragments 
to evaporate and mix with the bubble gas fairly rapidly, within say a 
few times 10* yrs. But for large galactic scale higher mass halos, I / 
is an order of magnitude larger, and tev is likely to be much larger. 

In order to examine if the effect of the RT instability signif- 
icantly changes the further evolution of the outflow, we have con- 
sidered the following simple toy model: We assume that the shell 
is broken periodically due to the RT instability during the accel- 
eration phase. At each such epoch we assume the shell looses a 
significant amount of its mass. However the hot bubble is still ex- 
panding outwards and sweeps up material from the surrounding to 
form a new shell. We assume the newly created shell has initially 
some small fraction (1%) of the original mass and expands with the 
same velocity it had before fragmentation. 

The effect of the Rayleigh-Taylor instability following the 
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for a halo of mass 10 M, 







above prescription is shown in Fig, 
collapsing at Zc — 4. This example would be relevant for Lyman 
break galaxies which are thought to be galactic scale objects having 
outflows. 

In panel (a) and panel (b) of Fig. [16] we compare the radius 
and velocity profile for the cases with and without the RT instabil- 
ity. The unperturbed solutions are shown by solid lines while the 
solutions taking account of RT instability as described above are 
shown by dashed and dotted-dashed lines. We have assumed that 
the shell gets fragmented periodically every 4 x 10^ yrs while it 
is accelerating. At such times the shell looses 99% of its mass and 
continues to sweep with the same speed as it had before fragmen- 
tation. We show two cases: one in which we add the fragmented 
cloud mass into the hot bubble assuming that they evaporate com- 



pletely (dotted-dashed lines) and a second in which the fragmented 
clouds do not add any mass into the hot bubble (dashed lines). 

We see from this figure (panel (b)) that the velocity in the ac- 
celerating stage are affected by the breaking of the shell due to the 
RT instability, however the final radius and velocity remain prac- 
tically unaltered. The maximum velocity achieved by the wind in- 
creases in the case where we take account of the RT instability, as 
the mass of the shell reduces periodically. The breaks seen in the 
velocity of the outflow mark the position when the shell fragments. 

The temperature of the hot bubble (panel (c) in Fig. [T6j also 
finally reaches the same asymptotic value. However, depending on 
whether the fragments can evaporate (dotted-lines) or not (dashed 
lines) the exact evolution of the temperature can be different. The 
fragment/hot bubble interfaces will have a range of temperatures 
from T ~ lO** — 10*^ K, with the metals from the hot fluid mixed 
in with the cooler shell gas. There could also be mixing due to 
secondary Kelvin-Helmoltz instabilities which operate at the inter- 
faces and/or the turbulence mentioned above. These interface re- 
gions could lead to a host of absorption lines in both the spectrum 
of the wind blowing galaxy and background quasars seen through 
the outflow region. 

From the above analysis we also see that the final radius and 
velocity of the outflows do not change significantly. Hence we 
will not consider the RT fragmentation while computing the global 
properties of the outflow. 



7 DETECTION OF METALS 

In the semi-analytic models described and explored in detail above, 
we can see that winds can efficiently transport metals into the IGM 
for a wide range of model parameters. The question we wish to 
ask now is that, what are the required physical conditions in this 
expelled gas that will enable us to detect them in spectroscopic 
observations. The detectability of the metals in a given ioniza- 
tion state depends on coUisional ionization, photoionization and 
recombination. Thus one needs to incorporate the effect of these 
processes into our calculations. The most easily detectable dou- 
blet transitions from the highly ionized gas are C IVAA1548, 1550, 
Si IVAA1393,1402 and O VlAA1032,1037. The respective ioniza- 
tion potentials are 64.5 eV, 45.1 eV and 138.2 eV. The ionization 
fraction for C IV, Si IV and O VI peaks at temperatures of 10"" K, 
6.4 X 10* K and 2.8 x 10^ K respectively for collisional ioniza- 
tion (Sutherland & Dopita 1993). As the recombination time-scales 
are larger for mean IGM density at z ~ 3 the gas temperature is 
not controlled by the ionization equilibrium and the wind dynamics 
plays an important role in the gas temperature and density. 



7.1 Ionization correction 

Now we calculate the ionization fraction for the three species 
namely C IV, Si IV and O VI as a function of density and tempera- 
ture of the medium. In order to compute the ionization corrections 
correctly we have used the spectral simulation code cloudy (ver- 
sion C06.02 - see Ferland et al. 1998 for details). We assume the 
gas to be a plane parallel slab with total hydrogen column density 
of 10^^ cm~^. We have used the QSO and Lyman break galaxy 
dominated meta-galactic UV background at z — 3 computed by 
Haardt and Madau (2001) as the ionizing source. The gas density 
is varied between 10~^ to 1 cm~^ and the temperature of the gas is 
varied between 10'*K to 10^ K. 
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Figure 16. The effect of RT instability on outflow. We show R^, Vg and of an out flow originating from a 10^^ Mq halo collapsed at 2^ = 4 in panel 
(a), (b) and (c) respectively. The solid curves are for unperturbed solution while the other two are the solution taking account of the RT instability. The dashed 
lines are obtained when we consider no evaporation of the fragments while the dotted-dashed lines are plotted when we assume complete evaporation of the 
fragments and instantaneous mixing of the fragmented mass into the hot bubble. 
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Table 1. Ionization fraction 
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Figure 17. Ionization fraction of C IV (top) and O VI (bottom) as a function 
of density of the medium. Different colour coded curves are for different 
temperature. 



The ionization fraction of C IV and O VI as a function of den- 
sity for a range of temperatures a are shown in Fig.[T7] In Table. [T] 
we give the fractional ionization of the three species (C IV, Si IV 
and O VI ) at three different temperature and densities. Here, p is 



Density 


species 


fractional ionization at T = 
10-* K IQS K 5 X 10^ K 


P 


C IV 
O VI 
Si IV 


1.6 X 10-2 
2.6 X 10-1 
5.6 X 10-^ 


4.8 X 10-1 
4.1 X 10-2 
3.6 X 10-^ 


1.0 X 10-5 

3.5 X 10-3 

1.6 X 10-* 


O.lp 


C IV 
O VI 
Si IV 


9.4 X 10-5 

1.5 X 10-2 
2.0 X 10-1" 


9.1 X 10-^ 
4.4 X 10~i 
1.0 X 10-12 


3.6 X 10-* 
2.1 X 10-5 
1.8 X 10-3 


lOp 


C IV 
O VI 
Si IV 


2.7 X 10-1 
1.7 X 10-1 
6.3 X 10-3 


3.8 X 10-2 
6.3 X 10-2 
7.2 X 10-^5 


3.3 X 10-* 
2.2 X 10-2 
9.7 X 10-'^ 



the mean IGM density at z — 3. We wish to point out that results 
will not depend too much on our choice of N(H) as long as N(H I) 
is optically thin at the Lyman limit. It is clear from the figure that at 
low densities and temperatures the effect of photoionization is very 
important. The ionization fraction of the two species are controlled 
by collisions (when it is independent of density) at higher densities. 
The vertical line in this figure marks mean IGM density at z = 3. 
We notice that the maximum Si IV fraction is 6.3 x IQ-^ when we 
consider T — 10'* K, p — Wp. For higher temperature and lower 
density gas the fraction of Si IV is negligibly small. 

From the table it is clear that when one considers the under- 
dense region (density ~ O.lp) ionization fraction of C IV is less 
than IQ-'' for T > 10^ K. This will make detection of C IV in ab- 
sorption towards high redshift bright objects virtually impossible. 
This exercise clearly states that if the metals are going to be dis- 
tributed in voids with average density less than the cosmological 
mean density then it will be very difficult to detect them in C IV 
absorption. The maximum ionization fraction for O VI we get for 
this case is 1.5%. The column density of O VI is given by. 



N(0 VI) = 7.4 X 10~* X N{H) 



fovi 



(21) 



0.01 J VO.Ol Zq^ 
Typically a low density region with N(H) = 102i cm-2 will produce 
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detectable O VI absorption when the gas is at IC* K. However, the 
detection becomes difficult if the temperature of the gas is higher 
than 10* K. Clearly when O VI is detected from the low density 
medium it will have signatures of photoionization and not that of 
coUisional ionization. 

If the over dense regions, like that probed by the high column 
density Lyman-Q forest, are polluted by metals then it will become 
easy to detect C IV and O VI absorption from this gas if the tem- 
peratures are close to photoionization temperatures. Again if the 
gas temperature is higher than lO^K it will become more difficult 
to detect C IV. Having got a rough idea of the temperature density 
ranges where O IV and C IV fraction is higher in the presence of 
meta-galactic UV background radiation, we now consider different 
cases in our models. 

Spectroscopy of bright QSOs, GRBs and Lyman break galax- 
ies allow us to probe the metals in the IGM through absorption line 
seen in their spectra. Direct spectroscopy of galaxies probes the 
ongoing winds in the galaxies where as the spectroscopy of QSOs 
allow us to probe the global enrichment of the IGM and enrich- 
ment around bright galaxies. Clearly different observations probe 
the nature of the wind at different stages of its evolution. 



7.2 Metals in free wind 

As seen from Fig.[T] the metal enriched gas after escaping the galac- 
tic disk will travel outwards in the form of free wind before encoun- 
tering the inner shock at _Ri . Metal absorption in the spectrum of 
galaxies with small radial separation could come from the free wind 
material coming out of the galaxy. In our model the free wind ma- 
terial extends upto the iimer shock {Ri) which can go beyond virial 
radius of the halo. When this gas is shocked at the inner shock the 
temperature goes up and hence detection probability of the metals 
goes down. Even a naive calculation (see below) then shows that 
the free wind material has sufficient column density for detection 
in line absorption. 

We assume that the free wind is not very hot. This may hap- 
pen because cooling is very efficient in the ejecta of individual SNe 
or due to subsequent adiabatic expansion of the thermally driven 
galactic wind. We integrate the density of the free wind from Ri to 
R\ to find out the column density of H, assuming for illustration a 
constant mass outflow rate. Since the free wind material is expand- 
ing asymptotically at a constant velocity the density will 
scale pn] {r) oc r~^. The normalization can be fixed at the initial ra- 
dius Ri from mass conservation to be pw{Ri) = / {A-kR^v^). 
Integrating this profile we get column density of hydrogen 



1 

Ri 



(22) 



where rrip is the proton mass. Taking Mn> ~ IOMq yr~^, Vm ~ 
400 km s~^, Ri = 1 kpc and assuming i?i >> Ri we get N(H) ~ 
2 X 10^" cm^^. From Eq. it is clear that O VI column density 
of this gas will be 1.5 x l(r^ cm~^ if the metallicity is 0.1 and 
fovi is 0.001. Note that in our models A/„ is a function of time 
and the above exercise gives us only an order of magnitude estimate 
of typical expected column density. 

In our models high mass outflow rates are expected in high 
mass galaxies as they have higher SFR. However, the exact value 
of /o vi (or fc iv) will depend on the metallicity and temperature 
of the free wind gas. Ideally we expect the initial temperature of 
the free wind to be equal to or less than the temperature of the su- 
perbubble that initiates the flow. This will be typically in the range 



few time 10^ K to 10*^ K (see also Fig 3. of Efstathiou 2000). Sim- 
ulations with single starburst in the disk suggest the temperatures 
could be as high as 10^ K in the inner side of the free wind (see 
Fig. 4 of Mac Low & Ferrara, 1999 and Fig. 6 of Fujita et al. 2004). 
If the free wind is at high temperature then detectability in the UV 
absorption is difficult. However if the temperatures are more like 
few 10^ K then we expect to detect the free wind in UV absorption 
lines discussed above. 

It is important to note that the outflows detected in the spectra 
of high redshift Lyman break galaxies, starburst and dwarf galax- 
ies at low redshifts also show standard low ionization interstellar 
absorption in the outflowing gas (Pettini et al. 2000; Rupke et al. 
2002; Martin 2005). Presence of Na I absorption in the outflow 
means, clearly some part of the outflowing gas is optically thick to 
the Lyman continuum photons originating either from galaxy itself 
or from the meta-galactic UV background radiation. Note that most 
of the observational papers favor the distance of the outflowing gas 
to be few kpc away from the galactic disks. If we assume that the 
outflowing material from the galaxy has properties very similar to 
that in ISM then we can envisage a situation in which the outflow- 
ing material has mixture of hot, warm and cold media (Heckman et 
al. 2000). In such a case a free wind material can produce absorp- 
tion lines covering wide range of ionization states (see also Murray 
et al. 2007). 

7.3 Metals in the bubble and shell 

Next we consider the bubble gas after the free wind enters the inner 
shock. As discussed before, the free wind material and the (1 — e) 
times the swept up material are in the bubble. The bubble density 
is usually low compared to the ambient medium. When the shell 
is well within the virial radius the bubble density is contributed 
by both the leftover gas from the swept up material (or evaporated 
material from the shell) and the shocked wind gas. In our models 
this is roughly in the range 10"^ to 10"^ cm . However, the gas 
temperature is usually higher than 10*' K. At these temperatures 
the gas ionization will be dominated by collisions and both C IV 
and O VI will not be the dominant ionization state of the respective 
atoms. 

It is clear that for a uniform bubble material it will be diffi- 
cult to detect the outflowing gas in the UV absorption. It is also 
important to realize that when the outflowing gas is near the virial 
radius SFR in the galaxy is also near its peak. Thus the study of ab- 
sorption lines in bright galaxies will correspond to the situation in 
which the shock front is either inside the virial radius of the galaxy 
or just about leaving it. Thus, if the UV absorption lines seen in 
the spectra of galaxies originate from the bubble material then we 
needs inhomogeneities in the bubble gas as well. This could ei- 
ther be from (i) the shocked wind material that already has clumps 
as we discussed above (ii) or the broken shell fragments from the 
RT instability. In the former case the high density ISM gas will al- 
ready have metals mixed into it. This will enable one to detect even 
low ionization material. However, in the latter case the broken shell 
material may predominantly have swept up material that may not 
be highly enriched. High ionization species can also originate from 
the conductive interface between hot metal rich bubble material and 
cold shell material (as discussed in the case of Weaver et al. 1997). 
From our discussion on RT instability in section 6, it is clear that 
such fragments survive for a longer time-scale against evaporation 
in high mass halos. 

When the shock front leaves the virial radius, the average bub- 
ble density is (1 — e) times the IGM density. For e = 0.9, then 
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typically, the bubble is underdense by a factor 10 compared to the 
uniform IGM. From our earlier discussion on individual halos we 
note that the bubble temperature in this stage of evolution is de- 
cided by adiabatic cooling and is some where in the range 
to 10^ K. From Table [T] it is clear that it will be difficult to detect 
both C IV and O IV in absorption leave alone detecting the low 
ionization lines seen in the outflow. If UV absorption is produced 
from the bubble gas we need inhomogeneities with higher density 
and lower temperature (as discussed above) compared to the mean 
bubble density and temperature. 

On the other hand the shell density is high enough and tem- 
perature is low enough (when the shock is well within the virial 
radius) so that the metals can be detected in absorption. The shell 
material is dominated by the swept up mass from halo gas/IGM. 
If the IGM is not pre-enriched by the metals then one would not 
detect the shell in metals line. 



7.4 Effect of pre-enrichment 

As discussed in the introduction, C IV, Si IV and O IV absorp- 
tion lines detected in the spectra of high redshift QSOs probe the 
metal enrichment history of the universe. The discussions presented 
above suggests that when the metals are expelled from the galax- 
ies they usually displace the material in the ambient medium into 
a thin shell and spread the metals into a low density high temper- 
ature bubble. In a simple uniform density model without inhomo- 
geneities it will be difficult to detect these shells through C IV and 
O VI absorption lines. 

Note that when an outflow is traversing through a medium 
which is already enriched by earlier generation of star formation 
the swept-up shell material can be detected in absorption lines. For 
example the column density of carbon at the virial radius of the 
halo is 



N{C) 




Here we have assume that the halo gas is pre-enriched with a metal- 
licity of 10^'^ Zq. When the shell is within the virial radius of the 
halo its density varies between 10^^ to 10"'^ cm~^ and the tem- 
perature is ~ 10* K. The C IV fraction will be close to 0.3 and the 
shell will produce detectable C IV absorption. However when the 
outflow is in IGM, the density of the shell material is very low and 
the temperature is more than 10^ K. This makes it difficult to detect 
the shell in C IV even though the column density of carbon is still 
of similar order compared to the case when the outflow is inside the 
halo. 

Interestingly our knowledge of metals in the Lyman-a clouds 
comes from observations of high density regions with over density 
in excess of 6 or so. One way of incorporating metal into these 
over dense regions without altering the statistical properties of the 
Lyman-Q forest distribution is to inject material into the IGM very 
early (while a. z = 3 over dense region was very close to the mean 
density at an earlier epoch) and allow the metal mixing through 
the evolution of gravitational perturbations that produce Lyman-a 
forest absorption lines. In this scenario it is quiet possible that the 
metal absorption lines we see in the high density Lyman-a forest 
without the signature of coUisional excitations could have come 
out of very early generation of low mass galaxies. Even though 
our models do not capture the non-linear evolution of the IGM, in 



the following section when we discuss the global properties of the 
outflows we will show that the IGM could be pre-enriched with 
metals. 



8 GLOBAL CONSEQUENCES OF OUTFLOWS 

From calculating the evolution of a suite of individual outflow mod- 
els, we can study several global properties of the wind affected re- 
gions. This is mainly to understand the effect of the outflows on the 
IGM. One simple quantity is the porosity Q{z), defined by adding 
up the outflow volumes around all the sources at any redshift: 

oo oo 

Here, d^N{M, z, z')/dz' dM comes from Eq. l|T]l and Rs comes 
from solving for the outflow dynamics. The lower limit A/iow is 
determined from the cooling criteria and the effects of reionization 
feedback, as described in section [2] For Q < 1 the porosity gives 
the probability that a randomly selected point in the universe at z 
lies within an outflow region. For Q ~ 1, it is more useful to define 
the associated filling factor of the outflow regions, which if out- 
flows are randomly distributed is given by, F = 1 — exp[— Q(z)]. 
These considerations ignore source clustering, which can be im- 
portant for rare sources and lead to a smaller F. Note that Q can 
exceed unity while F{z) < 1 always. 

We can also define porosity weighted averages of various 
physical quantities associated with the outflows, and their proba- 
bility distribution functions (PDFs). For any physical quantity say 
X, its porosity weighted average is defined as 

^^[i?s(l + ^)]'X (25) 

The cumulative PDF can be obtained by replacing X in Eq. l[25j by 
the Heaviside Theta function Q{X{M, z, z') — Xo) and the differ- 
ential PDF can be obtained by differentiating this respect to Xq. We 
compute below such porosity weighted averages and the associated 
PDFs of various physical properties related to the outflows. 

One of the main motivations to study the SNe driven galactic 
outflows is to understand the metal pollution of the IGM. To com- 
pute this, we calculate the mass of metals coming out of the galaxy 
through an outflow as described in Appendix |A] We divide the to- 
tal mass of ejected metals in a fiducial volume of the universe by 
the average baryonic mass in this volume and refer to this as the 
average global metallicity, Z of the IGM. 

Note that the global properties of the outflows will depend 
on the parameters describing star formation, reionization and cos- 
mology. As mentioned in section|2] we will consider both 'atomic 
cooling' and 'molecular cooling' models of Paper I, which are con- 
sistent with the observed high redshift UV luminosity functions of 
galaxies and observed constraints on reionization. 

8.1 Atomic cooling Model 

First, consider the model parameters for individual halos as dis- 
cussed in section [4] [i.e we take /* — 0.5, n = 1.0, fh = 0.1, 
e = 0.9, Ew = 0.1, 77 = 0.3, = 50 Mq with cosmological 
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Figure 18. The volume filling factor as a function of redshift for 'atomic 
cooling' model A. Total filling factor of the universe is shown by the solid 
line. The dotted-dashed line shows the contribution coming from the out- 
flows which are frozen into the bubble flow where as the dashed curve 
shows contribution by the rest of the outflows which are still moving at 
a supersonic velocities. 



parameters from the WMAP 3rd yr data]. We refer to this as our 
fiducial model A. 



8.1.1 The volume filling factor 

In Fig. [Is] we show the volume filling factor F{z) for model A, as 
a function of redshift. The solid line shows the net volume filling 
factor. We can also split this up into contributions from outflows 
which have already frozen into the Hubble flow, say Fh{z) and 
those which have not (and moving still at a supersonic velocity), 
denoted by Fa(z). The dash-dotted and dashed lines shows respec- 
tively Fh and Fa- At high redshifts (z > 6) the volume filling 
factor is dominated by Fa', i.e. the outflows which have not frozen 
into the bubble flow. However below redshifts of about 3, the con- 
tribution from the bubble frozen outflows starts to dominate. It is 
clear from the Fig.[T8]more than 30% of the universe is affected by 
the outflows even at z ~ 6 and this increases to ~ 60% by 2 ~ 2 
for these model parameters. 

It is also important to know the mass range of halos that con- 
tribute significantly to the volume filling factor at different epochs. 
Fig. [T9] gives the contribution to the volume filling factor by ha- 
los in different mass ranges. At redshifts z > 2 the filling factor 
is dominated by galaxies with mass range 10^ — 10^ Mq while 
higher mass halos start to dominate at lower redshifts. Earlier work 
by (Madau, Ferrara, & Rees 200 1^ in fact examined the effect of 
outflows from galaxies with a typical halo mass of W^h~^ Mq. 
We also find that such halos are important contributors to F. How- 
ever our work includes halos of all mass ranges that are allowed 
by our cooling criteria, and we find that this leads to a significant 
increase in the volume filling factor of outflows to give _F ~ 0.7 
by z ~ 1. It is also interesting to note that the galaxies contribut- 
ing to the observed high-z UV luminosity functions have typically 
masses M > IO^Mq (see Fig. 4 in Paper I). Therefore galaxies 




Figure 19. Contribution of halos with different mass ranges to the filling 
factor of IGM. Curves are drawn for mass range lO'^ — 10** Mq (dashed), 
10* - 10^ Mq (dot-dot-dashed), 10^ - 10^° Mq (dotted) and lO"^ - 
10^^ Mq (dash-dotted). The total filling factor is also shown by solid line. 



which dominate the volume filling factor are not yet detected di- 
rectly. 



8. 1.2 Porosity weighted averages and PDFs 

The porosity weighted averages of several physical properties of 
outflows are shown in Fig. [20] Panel (a) shows the volume filling 
factor, while panel (b) gives the porosity averaged dark matter mass 
of the halo from which the outflows are originating. Both these 
aspects of outflows have been discussed above. Panel (c) shows 
the porosity averaged comoving radius of the outflow in units of 
kpc. This increases from ~ 100 kpc at 2; ~ 10 to ~ 500 kpc at 
2 ~ 3 upto 1 Mpc by 2 = 1. More detailed information is available 
in panel (a) and (b) of Fig. |2T| where we show respectively, the 
differential and cumulative PDFs of the outflow radius (proper), for 

2 = 2, 3, 4 and 5. For clarity, we have shown the differential PDFs 
only for z — 5 and 3. It is clear from the figure that the median 
proper radius of the bubble is 60, 80, 110 and 200 kpc, respectively, 
for 2 = 5, 4, 3 and 2. These scales are perhaps just below the scales 
in the matter power spectrum probed by the Lyman a forest at 2 ~ 

3 (Croft et al. 1998). The outflows also fill a significant fraction 
of the volume, with F ^ 0.6 by 2 ~ 2. However the outflows 
are believed to propagate more into the voids than along filaments 
(Theuns et al. 2002), and so may not perturb the overdense regions 
from which the Lyman a forest absorption originates. We discuss 
this issue further below. 

In panel (d) of Fig. [20] we show the porosity weighted peculiar 
velocity of outflows. This velocity is ~ 75 km s~^ at high red- 
shifts 2 ~ 12, reflecting the fact that most of the outflows are very 
young at this stage. As time increases the average peculiar veloc- 
ity decreases to about 20 km s~^ at 2 ~ 3 — 4, as considerable 
volume is filled by the flows, originating from early generation of 
galaxies, that are either Hubble frozen or moving with much lower 
peculiar velocities. However, at 2 < 3 the porosity weighted pe- 
culiar velocity increases due to an increase in the contribution of 
outflows from high mass halos. Note that as the volume filling fac- 
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Figure 20. The global properties of outflows. Panel (a) shows the volume 
filling factor _F. The porosity weighted average dark matter mass (in Mq ), 
comoving radius (in kpc) and peculiar velocity (in km s^^) of the outflows 
are shown in panel (b), (c) and (d) respectively. The porosity weighted den- 
sity (in cm~^) and temperature (in K) of the hot bubble are shown in panel 
(e) and (f). In panel (e) we also show the mean baryonic density of IGM by 
dashed line. 



z = 4 z-2 





Rj ( kpc ) 



Rj ( kpc ) 



Figure 21. The differential (panel a) and cumulative (panel b) PDFs of the 
outflow radius (Rs) for model A. 



tor increases, more and more outflows will interact and the peculiar 
velocity of the outflows can then lead to supersonic turbulence in 
overlapping regions. Such turbulence can also lead to amplification 
of magnetic fields in the IGM, by the operation of the fluctuation 
dynamo (cf. Zeldovich, Ruzmaikin & Sokoloff 1990; Brandenburg 
& Subramanian 2005). 

The solid and dashed curves in panel (e) of Fig. |20] show re- 
spectively, the porosity weighted number density of the hot bubble 
and the mean baryonic density of the IGM. It is clear that the bub- 
ble density is always less than the mean IGM density at the same 
epoch by a factor ~ 10. This factor is determined by the value of 
the entrainment parameter, e, since the bubble mass is eventually 
dominated by the mass swept from the IGM and not by the mass 
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Figure 22. Probability distribution function (PDF) of temperature and den- 
sity of the hot bubble. In panel (a) we show the PDF of bubble temperatures. 
In panel (b) we show the cumulative distribution of bubble temperature. In 
panel (c) and (d) we show the similar distribution of bubble density. 



coming from the galaxy. As discussed in the previous section if the 
hot bubbles have low density then the metals in these bubbles will 
not be detected in the UV spectroscopy of the bright objects. As 
the fresh outflow most probably propagates into the voids the met- 
als detected in the high column density Lyman-a absorption lines 
have to come from the early generation of outflow that pre-enriched 
the IGM. 

Th e average bubble temperature is shown in panel (f) of 
It decreases from about 10^ K at z ~ 10 to ~ 10^ K 
3 — 4. More detailed information is available in panel (a) 



Fig. 
at z 



20 



and (b) of Fig. [22] where we show respectively, the differential and 
cumulative PDFs of the bubble temperature, for z — 2, 3, 4 and 5. 
For visual clarity, the differential PDF is given only for z = 3 and 
z — 5. At z — 5 , more than 70% of the bubbles are at a temper- 
ature higher than 10"" K where as by z = 3 this fraction decreases 
to less than 15%. This is mainly due to the adiabatic expansion 
of the bubbles. The void regions these bubbles fill will neverthe- 
less be at a higher temperature than the photoionized IGM. At even 
lower redshifts when halos with M > W^AIq start to contribute 
significantly to the volume filling factor, the fraction of bubbles at 
Tt > 10^ K increases, although more than 80% of the volume 
filled by the bubbles still have Tb < 3 x 10^ K (see the cumulative 
PDF). The photoionized IGM is expected to have temperatures in 
the range 2 — 4 x 10'* K in the redshift range 2 < z < A (Schaye 
et al. 2000). The model discussed above has 55% and 20% of the 
volume being filled by the gas with T > 5x 10'*Kat2 = 4 and 
3 respectively. Influence of such a gas to the observed properties 
of the Lyman-Q forest is negligible if the hot gas predominantly 
percolates into the low density voids. However, as pointed out by 
Theuns, Mo & Schaye (2001), if the gas around massive galaxies 
are uniformly heated to high temperatures their effect will be felt 
in the high column density end of the Lyman-a forest. 

In panels (c) and (d) of Fig.|22]we also show the correspond- 
ing PDFs of bubble metallicity. Most of the bubbles have metal- 
licities between 0.01 — 0.1 Zq, with the differential PDF peaked 
around the lower value. The average global metallicity evolution is 
also of interest. We show this in Fig. [23] for our fiducial model A. 
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Figure 23. The global average metallicity evolution of the IGM (solid) and 
galaxies (dashed) as a function of redshift. The Vertical dotted line corre- 
sponds to the redshift where Q = 1. 



The epoch when the porosity becomes unity is marked by a ver- 
tical dotted line. One can see that the globally averaged metallic- 
ity gradually builds up from about Z ~ 10^^ Zq^ at 2; ~ 10 to 
Z ~ 2 X 10~^ Zq at z ~ 2, by which time the porosity has just 
exceeded unity. Further, any halo collapsing after the epoch where 
(5 ~ 1, should have at least this amount of metals present in its ISM 
even before it starts forming stars. We will also compare below, the 
global metallicity evolution among different models. Fig. |23|also 
show the average of all metals produced in the galaxies. It is clear 
that for the model parameters we have chosen here considerable 
percent of the metals that are produced remains in the halos. Thus 
accurate estimation of metal budget at high-2: can be used to con- 
strain the models. It has been argued that only ~ 60% of the metals 
produced can be accounted for by the metals in the high redshift 
galaxies and > 40% of the global metals produced have to be ac- 
counted for by the metals spread outside the bright galaxies (see 
Bouche et al. (2007) for recent metal budget). This means we may 
need to expel more metals than what is the case with our fiducial 
model considered here. 



8.2 Reionization feedback 

In this paper we have calculated the global outflow properties tak- 
ing into account self-consistently the feed back due to the reion- 
ization history. This is an important improvement as most of the 
earlier works on galactic outflows do not take into account the ra- 
diative feedback in the star formation models. In Fig. [24] we high- 
light this effect by comparing the result of our self-consistent model 
with two toy models with sudden reionization. The solid line cor- 
responds to the evolution of F{z) in Model A calculated with the 
self-consistent ionization history as in Paper I (See their Table 2 
and Fig. 2). In this model the reionization occurs at z^e = 7 and 
the universe gets 50% ionized at 2: ~ 8.5. 

The dashed line assumes that the universe underwent an abrupt 
reionization at z^e = 6, while the dash-dotted line assumes Zre — 
11. We see that early abrupt reionization at Zre = 11 leads to a 
significant fall in the volume filling factor, especially at z > 2. This 
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Figure 24. Effect of reionization on the volume filling factor. The volume 
filling factor for model A is shown by solid line. The dashed and dotted- 
dashed lines correspond to a step reionization model with Zre = 6 and 
Zre = 11 respectively keeping all other model parameters identical to that 
of model A. 



is because the smaller mass halos which contribute dominantly to 
volume filling the IGM are suppressed from forming stars due to 
early reionization. Below this redshift, higher mass halos start to 
contribute leading to a rise in F{z). On the other hand having an 
abrupt reionization at a latter redshift Zre = 6, leads to 90% of 
the volume being filled with outflow material by z = 3. Therefore 
a determination of the volume filling factor due to outflows is a 
sensitive probe of the ionization history. This point has also been 
independently made recently by Fieri et al. (2007). 



8.3 Global properties for a range of model parameters 

In Table |2] we have listed a few models where we vary some pa- 
rameters around that adopted in the fiducial model A. We wish to 
investigate the effect of this change on the global outflow prop- 
erties. These models are all consistent with the observed high-z 
UV luminosity functions below z = 6. We also have taken the 
self-consistent reionization feedback as in Paper I. In Fig. |25] we 
show the global properties of the outflows for models listed in Ta- 
ble |2] In this figure we only concentrate on porosity (Q) (bottom 
panels), porosity weighted average temperature of the hot bubble 
(middle panels) and the global metallicity of the IGM (top panels). 
The other physical quantities, which are shown in Fig. [20] do not 
vary significantly while changing the model parameters. 

Model B examines the effect of adopting a lower /, = 0.25 
compared to model A. In order to still fit the UV luminosity func- 
tions of high redshift Lyman break galaxies, one then needs to take 
a lower value of k = 0.5 (Paper I). The porosity weighted mean 
temperature of the bubble for model B is slightly lower than that of 
model A. However, both the porosity and the average global metal- 
licity Z of the IGM are smaller than in model A. The porosity in 
fact never reaches unity in this case. The smaller values of Z and 
Q are the manifestation of the simple fact that total star formation 
in model B is half of model A. As discussed before the parameters 
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Table 2. Model parameters for our different models. 



model 






fesc 


OS 


71s 


V 




Z(Zo) 


Remarks 


A 


0.50 


1.00 


0.10 


0.75 


0.95 


0.30 


0.10 


0.01 


standard 


B 


0.25 


0.50 


0.10 


0.75 


0.95 


0.30 


0.10 


0.01 


change in /, & k wrt A 


C 


0.25 


0.50 


0.10 


0.75 


0.95 


2.00 


0.10 


0.01 


change in rj wrt B 


D 


0.50 


1.00 


0.10 


0.75 


0.95 


0.30 


0.02 


0.01 


change in wrt A 


E 


0.25 


1.00 


0.10 


0.85 


1.00 


0.30 


0.10 


0.01 


change in /* & cosmology wrt A 


F 


0.50 


1.00 


0.10 


0.75 


0.95 


0.30 


0.10 


0.03 


change in Z wrt A 




of this SFR model may correspond to a model with strong negative 
feedback from the SNe. 

In model C we explore the effect of increasing the the mass 
loading factor to 77 — 2. Note that we can not do this in model A 
because it assumes /* — 0.5, and such a large rj would require 
more baryons to be expelled than the galaxy originally contained. 
Hence in model C we assume all parameters as in model B, except 
for rj. Two major differences arise due to a larger rj now, compared 
to model A or model B. The average temperature of the bubble is 
significantly lower, and also the average global metallicity is higher. 
In section 5.2 (see also Fig.|10|), we have shown that increase in 77 
reduces the bubble temperature, prevents the escape of winds from 
massive halos and increases the asymptotic metallicity. The temper- 
ature is lower because in the initial stage of the outflow evolution. 



an increased mass loading (by a factor 2/0.3), leads to a signifi- 
cantly reduced temperature, and this in turn reduces the porosity 
averaged temperature. The increase in average global metallicity is 
also due to the increase in rj. From Eq. l |A3| l one can compute the 
mass of metals coming out of the galaxy. For 77 = 2, it is factor of 8 
higher than for 77 = 0.3, assuming the same /» = 0.25. The poros- 
ity curve does not show an increase al z < 2 (as seen in model 
B) due to the fact that the winds from the high mass halos do not 
escape efficiently due to higher cooling rate in the halos. 

We explore the possibility of having a lower wind efficiency of 
En, — 0.02 in model D. All other parameters are as in model A. In 
section 5.3 (see Fig.[TTJ we have shown the decreasing the energy 
efficiency (i.e changing either of IMF or e^) reduces the radius 
outflow radius. If the change in the efficiency factor is more than 
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2.5 we see that winds will not escape from the high mass halos. As 
the adopted change in for model A is factor 5 of that used in 
D we do not expect all the halos that have outflows in A to have 
one in D as well. As expected both the porosity and the bubble 
temperature are lowered. On the other hand it is interesting to note 
that the average global metallicity does not change for 2 > 6. This 
is because every halo which had an outflow in model A still has one 
in model D at these high redshifts. And the mass of metals ejected 
remains the same with redshift. However, at z < 6 the outflows do 
not escape from some higher mass halos due to lower and hence 
reducing the mass of ejected metals at these epochs. In any case it 
is interesting to note that the models with e„ < 0.02 do not achieve 

g = 1. 

In model E we change the cosmological parameters i.e. erg 
and Us. We take erg = 0.85 and = 1.0, more in tune with the 
WMAP 1st year data. We self-consistently calculate the SFR and 
reionization for this model. To fit the observed luminosity func- 
tion we need /, = 0.25 with n = 1. The reionization occurs at 
Zre = 8.4 with an electron optical depth of Te ^ 0.096. In this 
model, the porosity is initially (at z; > 3.5) higher than for model 
A, primarily because of the larger abundance of halos due to the 
higher erg and n^. But at lower redshifts the porosity decreases 
below that in model A, perhaps due to the lower The poros- 
ity weighted temperature and metallicity are close to that seen in 
model B confirming the above reasoning. 

We have also explored the effect of changing the bubble metal- 
licity in calculating the cooling of the bubble gas (model F). This 
has negligible effect on the global properties of the outflows. 

The most interesting outcome of the exercise presented in this 
section is that Q > 1 only for our fiducial model. Thus if star for- 
mation activities are sustained only in the atomic cooled halos then 
the parameters of the models should be close to that of our fidu- 
cial model in order to completely fill the IGM with metals. Overall 
the alternative models lead to a lower porosity and volume filling 
factor, as well as lower bubble temperatures. The average global 
metallicity can be larger than in model A if the mass loading factor 
rj is higher, but for the other alternatives Z is generally lower. 

Presence of radiative feedback makes it difficult to predict the 
trend of various quantities discussed above when we add molecular 
cooled halos in our models. This needs to be explored in a self- 
consistent way. This is what we do in the following sub-section. 

8.4 Models including Molecular Cooled halos 

In Paper I we considered a number of models where the lower mass 
cut-off of a halo which can host star formation is decided by assum- 
ing efficient molecular cooling (Tvir > 300 K). Such small mass 
halos are not detectable directly. But their influence can be felt indi- 
rectly via their ionizing efficiency. Several of the molecular cooling 
models of Paper I are also consistent with the available constraints 
on reionization. Outflows from such molecular cooled halos could 
also affect the properties of the IGM in important ways. Here we 
study the outflow properties of two such models, namely Model 
M2 and M3 ( see Table 3 of Paper I). In model M2 we assume that 
a fraction /, = 0.1 of the baryons in a molecular cooled halo is 
turned into stars having a normal salpeter IMF. This is in addition 
to the star formation in the atomic cooled halos, where the parame- 
ters for star formation are as in model A. Such a model not only fits 
the high-z UV luminosity function, but also has an electron scat- 
tering optical depth ~ 0.105, consistent with WMAP 3rd year 
data. The model leads to a complex ionization history, with the final 
reionization at Zr^ = 5.9 (Fig. 2 of Paper I). 



1 F 




Figure 26. The volume filling factor as a function of redshift for molecular 
cooled model (model M2 of Paper I). The solid line shows the total filling 
factor of the universe. The dashed curve is the filling factor by the outflows 
which have frozen into the hubble flow of the IGM where as the dotted 
curve is the filling fraction by the rest of the outflows. 



In Fig.|26]we have shown the predicted volume filling factor 
of outflows obtained in this model. Comparing it with Fig. [TS] of 
model A, we see that star formation in the molecular cooled halos 
leads to significant volume filling of the universe, with a porosity 
Q{z) — 1, at even at 2; ~ 8. In fact, the hubble frozen outflows 
dominate below this redshift and the Q contributed by such out- 
flows itself becomes of order unity below 2 ~ 6. Subsequent struc- 
ture formation below z ~ 6 can lead to a significant fraction of 
the metal enriched outflowing material from the molecular cooled 
halos being incorporated into the mildly overdense regions of the 
IGM (relevant for the Lyman-a forest lines), and lead to a metal- 
licity floor of the IGM. 

' we show Z evolution in the molecular cooled model 
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(solid line) and compare it with atomic cooled model (model A: 
dashed line). It is evident that at any particular epoch the average 
global metallicity produced by these models do not differ much. 
However, as pointed out above, the porosity of outflows is larger, 
at an earlier redshift, in the molecular cooled model. In this model, 
any non-linear or mildly non-linear structures that are formed after 
z ^ 8 will have an initial metallicity of Z > 10^"* Zq even before 
the new stars are formed. Thus a normal mode of star formation in 
such small mass halos is a very effective means of spreading metals 
in to the mildly non-linear regimes at low redshifts that are traced 
by the Lyman-Q forest. Note that one can increase the metallicity 
of the IGM at 2 ~ 8 by changing our model parameters. For exam- 
ple, we get Z ~ 10^'^'^Zq when we use t] — 2 and /* — 0.1 in 
the molecular cooled halos. Note that this is very close to the mini- 
mum metallicity suggested by observations of Songaila (2001). In- 
creasing /, to 0.2 will increase the metallicity in the outflow by 
a factor of 4.6. However, if we wish to preserve the ionization his- 
tory then the UV escape fraction /e should be lowered to keep /, 
conserved. Note the metallicity of the IGM can be increased if we 
relax the condition of uniform mixing (in Appendix |A| and allow 
the wind to have higher metallicity than the average ISM. 

To compare the predictions of this molecular cooling model 
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Figure 27. Global average IGM metallicity as a function of redshift. The 
solid line is for molecular cooled model and dashed line is for atomic cool 
model, model A. The points indicate when the porosity, Q, becomes unity. 
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28. The global properties of outflows for molecular cooling model 
Paper 1). 



(M2) with the atomic cooling model A, in more detail, we show 
the corresponding porosity averaged physical characteristics of the 
outflows in Fig. [28] The panels in this figure correspond to the same 
quantities as in Fig. [20] Comparing Fig. [28] with Fig.|20] we see 
that in the case of the molecular cooling model M2, on the aver- 
age, the contribution to the filling factor comes from very low mass 
halos, with mass ~ lO'^ Mq at z = 8 to ~ 10^ A'/q at z = 3. 
For such low mass halos, the outflow radius is small on average. 




Rs (kpc) 

Figure 29. Porosity contributed by outflows of different sizes. The thick 
line shows the distribution for atomic cooling model A where as the thin 
line is for molecular cooled model M2 . 



But their number density is high enough that they can significantly 
pollute the IGM. The porosity averaged comoving size of the out- 
flows ~ 300 kpc at 2 = 3, is somewhat smaller than for model 
A. This can also be seen from Fig.|29|where the distribution of the 
bubble radius in molecular cooled model M2 is compared with the 
atomic cooled model A. The suppression of the large size bubbles 
(~ 100 kpc) in M2 is due to a larger radiative feedback than in 
model A. We also note that the median bubble size in the model 
is < 40 kpc (proper) for z > 3. This is also a factor two lower 
than that of model A. The average density of the hot bubble does 
not show any significant difference between atomic and molecu- 
lar cooling models as the density of the bubble is mainly govern 
by entrainment parameter, e. One of the major differences is the 
lower porosity averaged bubble temperature. Since star formation 
and hence the total number of SNe is less in smaller mass halos the 
hot bubble ends up with a lower temperature. Also the outflows pro- 
duced by these smaller halos lead to smaller peak velocities. And 
the resulting porosity averaged peculiar velocity for the molecular 
cooled model is much smaller at z > 3, compared to that of atomic 
cooling models. Low values of porosity weighted temperature, ra- 
dius and the peculiar velocity coupled with Q > 1 for the bubble 
makes this model more favorable to pollute the IGM without dis- 
turbing the observable properties of the Lyman-a forest. 

Many earlier works related to molecular cooled halos, assume 
a top heavy mode of primordial star formation in molecular cooled 
halos. For a top heavy IMF one may expect high mass SNe to dom- 
inate increasing the mass of ejected metals. It would be interesting 
to examine if these features can lead to a larger metallicity floor 
or a larger volume filling factor. To examine the effect of such a 
top heavy mode, we consider the model M3 in Paper I. This model 
assumes that molecular cooled halos form stars in the mass range 
50 — 500 Mq in a salpeter IMF. In this model the reionization 
occurs at Zre = 1 1.6 with — 0.155 slightly higher than that con- 
strained by the WMAP 3rd year data. 

We follow Furlanetto and Loeb (2003) and assume that in this 
case, one SNe explodes for every 460 Mq of star formed, with an 
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Redshift Redshift 

Figure 30. The volume filling factor and global average metallicity pre- 
dicted by 'top heavy' molecular cooled model. Panel (a) shows the volume 
filling factor where as panel (b) shows the average global metallicity of the 
IGM. In panel (a) we also show the break up Fj\ (dashed line) and Fh 
(dotted-dashed line). 

energy output of lO^'^ ergs, and ejecting 4 Mq of carbon. We show 
in panel (a) and panel (b) of Fig. |30] the evolution of the volume 
filling factor and average global metallicity respectively, for this 
model. It is clear that even at high redshifts z ~ 10, outflows fill a 
significant fraction F ~ 0.25 of the volume. However the volume 
filling factor hovers around this value even till z ~ 3, and only after 
this epoch F continues to increase. This is because the universe is 
reionized at a very early epoch {zre ~ 11.6) in this model. Subse- 
quently, reionization feedback leads to a strong suppression of star 
formation in dwarf galaxies below this high redshift. 

This suppression also affects the average global metallicity. 
Even though each SNe in this model is producing 40 times more 
carbon than for model A/2, for a given amount of star formation 
the excess metallicity in M3 only by a factor 4. However, M3 pro- 
duces 17 time more UV ionizing photos compared to that of M2. At 
higher redshifts M3 has slightly higher metallicity. However at low 
redshifts higher mass unsuppressed 'atomic cooled' halos domi- 
nate in the global properties and hence the predicted average IGM 
metallicity for this model is similar to other models discussed ear- 
lier. The volume filling factor at low redshifts is moderately smaller 
than in model M2. 

Hence from above results we can conclude that the inclusion 
of star formation in molecular cooled halos will help in enriching 
the IGM at higher redshifts, and lead to a metallicity floor, pro- 
vided the molecular cooled halos are less efficient in reionizing the 
universe (for example, if stars form still with a normal IMF or the 
escape fraction of the UV photons are very low). 



9 DISCUSSION AND CONCLUSIONS 

We have examined in some detail here semi-analytic models 
of galactic outflows and their consequences for the intergalactic 
medium. Our models are constrained by available observations of 
both star formation and reionization. The nature of star formation 
is one of the key elements which controls the energetics of galactic 
outflows and also their metal content. At the same time, it also de- 
cides the reionization history which is an important input for the ra- 
diative feedback that suppresses star formation in low mass halos. 
We improve on earlier semi-analytical modeling of galactic out- 
flows in several important ways. Galaxies form stars continuously 
with the duration and fraction of baryons going into stars being 
constrained by the observed high-z UV luminosity functions (see 
Paper I for details). We take into account of existing constraints 
on reionization and most importantly, the implied radiative feed 
back on the star formation in a self-consistent manner We adopt 



the WMAP 3rd year cosmological parameters. We model galac- 
tic outflows in a manner similar to stellar wind blown bubbles (cf. 
Weaver et al. 1977), following the dynamics of both the outer shock 
Rs and also a possible inner reverse shock at Ri . This can naturally 
incorporate possible smooth transitions from pressure to momen- 
tum driven outflows. We use the modified PS formalism of ( |Sasaki| 
|1994[ l formalism to calculate the formation rate of dark matter ha- 
los and the global consequences of outflows, instead of taking just 
the time derivative of the PS function, which does not account for 
the destruction rate of halos. 

Outflows generically accelerate initially due to the increasing 
energy input from the galaxy and decreasing halo density profile. 
The outer shock velocity reaches a peak value which increases with 
halo mass, typically ranging from ~ 100 — 400 km s~^ for halo 
masses ~ 10^ — lO^^Af0 respectively. This phase lasts for a dy- 
namical time-scale after which the outflow decelerates, till it be- 
comes subsonic and freezes to the Hubble flow. The hot bubble of 
shocked wind material, has initial temperatures > lO*" K but sub- 
sequently cools due to adiabatic expansion. If there is significant 
mass loading from the galaxy, it can also cool radiatively to transit 
to a momentum driven flow. However, for most model parameters 
here, this does not occur The swept up shell gas typically cools 
efficiently while the outflow traverses the halo. During the acceler- 
ation phase the thin shell is also unstable to R-T instability leading 
to shell fragmentation (see also Ferrara & Ricotti (2006)) without 
significantly altering the final outflow radius. Such shell fragmen- 
tation could enhance the detectability of metal lines from outflows 
by providing evaporating interfaces where the hot metal-enriched 
bubble gas mixes with cool, dense shell material. 

During its evolution, individual outflows can travel well be- 
yond the virial radius of the host halo to proper distances ~ 200 — 
1000 kpc for the above halo mass range. The inner shock at Ri 
keeps up with the outer shock while the galaxy is actively form- 
ing stars, with typical values of R\/ Rs ^ 0.4 — 0.6. This is very 
similar to the wind structure seen in the simulations of Fujita et 
al. (2004). By the time the outflow becomes subsonic and freezes 
onto the Hubble flow, both the hot bubble and shell temperatures 
are ~ 10* K (determined now by photoheating), and the shell is 
likely to fragment and mix with the bubble and IGM gas. 

We explored in some detail the dependence of outflow prop- 
erties on the assumed initial conditions and various model parame- 
ters. We show that the initial conditions play very little role in de- 
ciding the nature of the outflows. We also check this by comparing 
our model predictions with scale-free solutions. 

We find that outflows can generically escape from the low 
mass halos {M < IO^A/q), that dominantly contribute to the vol- 
ume filling of the IGM. For galactic scale halos, as expected, having 
higher halo density /h, higher wind mass loading 77 or lower energy 
input efficiency (e^i/f,) makes it more difficult for outflows to es- 
cape. A burst mode of star formation generically leads to a smaller 
outflow radius, even for the same values of other parameters. As 
the outflow properties of low mass halos are less sensitive to above 
mentioned parameters the nature and efficiency of star formation in 
these objects decide the feedback due to galactic winds. 

The detection of metals from the outflows and in the IGM is 
one of the crucial issues in our paper. The gas phase metallicity in 
the ISM and the wind are self-consistently computed for a given 
star formation rate and IMF assuming instantaneous uniform mix- 
ing in the ISM. The detection of the expelled metals either in dif- 
ferent stages of the outflow or in the IGM depends crucially on the 
ionization state of the gas. Using photoionization calculations per- 
formed using Cloudy we have shown that the metals in free wind 
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and low density bubble will be very difficult to detect through stan- 
dard UV absorption lines of C IV and Si IV and O VI. We need 
clumped high density gas either coming from the ISM in the form 
of free wind or from the R-T instabilities to detect these absorp- 
tion lines. In particular some of these high ionization species can 
be detected in the conductive interfaces between the cold clumps 
and hot bubble material. We show that the metals in the underdense 
regions (p ~ O.lp) will be very difficult to detect if the tempera- 
tures are higher than 10''' K. Thus fresh outflow from galaxies that 
enter the IGM will be very difficult to detect. C IV and O VI ab- 
sorption lines are easily detectable if they originate from overdense 
regions (p > IQp) that are typically probed by high column den- 
sity Lyman-a absorption line with T ~ 10'* K. This is possible if 
IGM is already filled with pre-enriched gas by the time these over 
densities were order unity fluctuations (say z > 8). 

One of the important issues is to understand is how outflows 
impact on the physical properties of the intergalactic medium. We 
therefore computed the volume of the IGM affected by outflows, 
the porosity weighted averages and PDFs of several important out- 
flow characteristics, for a number of atomic and molecular cooling 
models. 

For our fiducial atomic cooling model A, more than 30% 
of the universe is affected by the outflows even at z ~ 6 and 
this increases to ~ 60% by 2 ~ 2. Galaxies with mass range 
lO'' — 10® Mq dominantly contribute to the volume filling fac- 
tor; higher masses only dominate at low redshifts. This is consis- 
tent with the suggestion of Madau, Ferrara and Rees (2001), of the 
dominant influence of halos of ~ IO^A^q in filling the universe 
with outflows; our work however includes halos of all mass ranges. 
Further, these galaxies which dominate in filling the IGM are not 
yet detected directly in the high-z UV luminosity functions (see 
Fig. 4. in Paper I), The porosity averaged outflow comoving radius, 
peculiar velocity and bubble temperature, evolve from ~ 100 kpc, 
~ 60 km s~* and 10^ K, respectively, at 2 ~ 10 to ~ 500 kpc 
20 km s~^ and ~ 10^ K at z ~ 3. The median value of the 
outflow radius is 60 kpc (proper) at 2; = 5 which increases to 
100 kpc at 2 = 3. Whereas at 2: = 5 , more than 60% of the 
volume filled by bubbles are at a temperature higher than 10^ K 
this fraction decreases to less than 15% by z = 3. This is mainly 
due to the adiabatic expansion of the bubbles. The void regions 
these bubbles fill will nevertheless be at a higher temperature than 
the photoionized IGM. Most of the bubbles have metallicities be- 
tween 0.01 — 0.1 Zq, with the differential PDF peaked around the 
lower value. The average global metallicity evolution is also of in- 
terest. For our fiducial model, this gradually builds up from about 
Z ~ 10"^ ^0 at z ~ 10 to Z ~ 2 X 10"^ at 2 ~ 2, by which 
time the porosity has just exceeded unity. 

We have examined several other atomic cooling models which 
are all consistent with the constraints on star formation obtained in 
paper I. These models also lead to significant filling of the IGM at 
z ~ 3 with metals (with -2.5 > [Z/Zq] > -3.7), the actual 
extent depending on the efficiency of winds, the initial mass func- 
tion (IMF) and the fractional mass that goes through star formation 
and cosmological parameters. The reionization history has a signif- 
icant effect on the volume filling factor, due to radiative feedback 
(see also Pieri et al. 2007). Further, a large fraction of outflows at 
2 ~ 3 are supersonic, hot (T > 10^ K) and have low density, mak- 
ing metal lines difficult to detect. These models may also result in 
significant perturbations in the IGM gas on scales probed by the 
Lyman-Q forest. 

On the other hand, we find that models including star forma- 
tion in molecular cooled halos with a normal mode of star forma- 



tion (or a lower UV escape fraction) can potentially volume fill the 
universe at z > 8 without drastic dynamic effects on the IGM, 
thereby setting up a possible metallicity floor (—4.0 < [Z/Zq] < 
—3.6). In fact, the bubble frozen outflows dominate below this red- 
shift and the Q contributed by such outflows itself becomes of or- 
der unity below z ~ 6. On the average, very low mass halos, with 
mass ~ 10*^ - 10^ Mq at z = 8 to ~ 10* - lO'' at 2 = 3 
dominantly contribute to the volume filling factor. The porosity av- 
eraged comoving radius of the outflows is less than 100 kpc at red- 
shifts where the Q ~ 1. The bubble temperature and peculiar ve- 
locities are also smaller than for the fiducial atomic cooling model. 
The above features make this model ideal to spread metals into 
the regions which will subsequently collapse to form the Lyman 
alpha forest regions, without unduly disturbing these regions dy- 
namically. To some extent this scenario is the extrapolation of the 
scenario of Madau, Ferrara and Rees (2001) which was applied to 
outflows from dwarf galaxies, to even lower mass halos. Interest- 
ingly, molecular cooled halos with a "top-heavy" mode of star for- 
mation are not very successful in establishing the metallicity floor 
because of the additional radiative feedback, that they induce. 

As we discussed above we use a functional form for SFR 
(constrained by the observations of high-z UV luminosity func- 
tions) without doing self-consistent calculations. However, to get 
reliable results from the self-consistent calculations one needs to 
be very clear about various physical processes that are involved. 
This is reflected in the fact that two recent simulations addressing 
this issue conclude differently. Scannapieco et al. (2006) find that 
supernova feedback and the resulting outflows decrease the frac- 
tion of baryons which are turned into stars, by factor ranging from 
2-4 as one changes the mass of the object, and also a related more 
rapid fall in the star formation rate with time. However, Koboy- 
ashi et al. (2007) find somewhat different results. They note that 
the two effects of supemovae, the increased metal line cooling due 
to the chemical feedback and the increased heating due to the en- 
ergy feedback both have opposite effects on the star formation rate 
in the galaxy. And indeed these two effects seem to cancel to pro- 
duce no net effect due to supemovae feedback on the star formation 
rate (see their Fig. 2 and the discussion). Therefore it is not entirely 
clear from these works the extent to which one needs to change /* 
and K. If SNe produce strong feedback effects as found by Scan- 
napieco et al. (2006) the star formation rate will fall sharply as a 
function of time. This corresponds to our models with low k that 
requires low value of /, in order to reproduce the UV luminosity 
functions. 

Future improvement of our work would involve replacing 
Eq. ^ with a model of star formation in a multiphase ISM in- 
cluding various heating and cooling processes, possible effects of 
dark halo clustering on the outflow properties and importantly set- 
ting up our semi-analytical model in the framework of a large-scale 
structure simulation. 
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APPENDIX A: METALLICITY EVOLUTION 

We calculate the metallicity of the bubble material as follows (also 
see Binney & Tremaine 1994). Suppose at any given time 5Ms is 
the amount of mass goes into star formation. Let's assume at that 
instance the metallicity of the IGM is Z. The amount of gas mass 
lost from the ISM is 5Mg = —(1 + i])SMs where we take the 
mass loss rate due to wind is rjSMs and we neglect the time delay 
between the star formation and SNe explosion. If one takes that p 
is the amount of heavy metals ejected per solar mass of star formed 
then the change of metal mass in the ISM is given by SMh — [p~ 
Z{l + ri)] SMs . On the other hand the increase of metals in the wind 
material is 5mh ~ ZrjSMs. Now Z = Mh/Mg. Differentiating 
this and substituting for 5Mg and SMh from above, one can get 
5Z = {p/Mg)SMs. Taking the initial gas mass in the ISM as Mo, 
we get dZ/dMs = p/(Mo - (1 + »7)Ms). Integrating this with the 
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boundary condition that Mh = when AIs 
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Using this relation we obtain the metal mass in the wind material is 
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In the asymptotic limit when all the star formation is over i.e. 

Ms /Mo = /. we get 



rrih = 



VP 



;Mo[il + v)f* 
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+ (l-(l + r,)/.)ln(l-(l + r,)/.)] 
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To give a rough idea of numbers involved, we take 77 — 0.3, 
/. = 0.5 and p = 0.1/50 (one SNe will form per 50 Mq of star 
formation and 0.1 Mq carbon will produce from each SNe). In the 
asymptotic limit this will give the metallicity of the wind material 
as 0.2Zq. For t] = 1.0 and /* = 0.25 the metallicity is 0.1^0. 
However, in realistic situation this gas is also going to mix with 
some fraction of halo/IGM material reducing the metallicity of the 
hot bubble. Another thing to note is that if the ISM of the galaxy 
is already enriched with metals then the amount of heavy elements 
transfered to the hot bubble from the wind is 



rrih = 



VP 



;Mo 



(1 + 77)2 

In fl- (1 + 77) 



(1+77) 

Af 



Mo 



Ms^ 
Mo 

+ ZovM, 



1-(1 + ^) 



Ms^ 
Mo 



(A4) 



3.5 



2.5 - 



2 - 



Density 
Velocity 
Pressure 



\ 



1.5 



0.85 



0.9 



0.95 



r/R, (t) 



Figure Bl. The shell structure for our model parameters when outflow is 
traversing inside the halo. We show the dimensionless scaled density, ve- 
locity and pressure as solid, dashed and dash-dotted line respectively. We 
show the profile upto the contact discontinuity. 



material and which in turn will enhance the density further, making 
the shell much thinner. 



where Zo is the metallicity of the ISM when the galaxy is formed. 



APPENDIX B: STRUCTURE OF THE SHELL IN 
ADIABATIC REGIME 

We consider here the inner structure of the swept up shell. This is 
important for detection of the shell material and for determining its 
cooling efficiency. Inside the halo, the density falls approximately 
as pB oc r^^ ® and L{t) oc t which leads to Rs{t) oc t^^^'^ . In 
this limit one can obtain a self-similar solution for the shell struc- 
ture following Weaver et al (1977) and Koo & McKee (1992). We 
solve the continuity equation, momentum conservation equation 
and energy conservation equation assuming a self-similar solution 
(c.f. Eq. (2)-(8) in Weaver et al. (1977) and Appendix B of Koo 
& McKee (1992)). The boundary conditions are obtained from the 
shock jump conditions assuming a strong shock. Fig. |Bl| gives the 
self similar scaled structure of the shell. We show as solid, dashed 
and dash-dotted lines, respectively, the dimensionless scaled den- 
sity p, scaled velocity v and scaled pressure P as a function of 
similarity variable A = r/Rs. The dimensionless scaled param- 
eters are defined as p(r) — p\{t)p{X), v{r) — vi{t)v{X) and 
P(r) = Pi{t)P{X) where the pi(t), vi{t) and Pi{t) are the post 
shock values given by the usual Rankine-Hugoniot jump condi- 
tions. We show the density, velocity and pressure upto the contact 
discontinuity, whose location is obtained from the criteria given in 
Koo & McKee (1992). From the figure it is clear that when the 
outflow is traversing through the halo the shell can indeed be ap- 
proximated as a thin shell. It is only ~ 10% of the radius. Within 
the shell the density can be very high. This help in cooling the shell 
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